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ABSTRACT 


Effects or rural land-use change on micro-climate 
in Central Alberta are investigated by comparing the energy 
‘balances of six differently vegetated small plots. The plots 
are located at the University of Alberta Agricultural Re- 
search Station at Ellerslie, 10 km south of Edmonton. 

The six surfaces are trembling aspen (Populus trem- 
uloides Michx.), white spruce (Picea glauca Voss), western 
snowberry (S. occidentalis), barley, fallow, and short grass. 
All plant species are common locally. The plots are situated 
in close proximity to one another, on a relatively flat plain, 
thereby negating the effects of varied topographic exposure, 
soil type, and timing of weather events on energy balance 
differences calculated among them. Any significant variation 
inegelmarads aivlOnsmecOlteheatetl is misensyDLemneat a1 lux, Or 
latent heat flux, can be traced to the influences of the 
landscape and vegetation surfaces themselves. 

Statistical analysis of the results reveals signifi- 
cant differences in net radiation, sensible heat flux, and 
soil heat flux. There is no significant difference in 
latent heat flux between the plots, despite wide variations 
in soil moisture content. Poplar has the highest net radi- 
ation values, while barley and grass record the lowest. 

The remaining surfaces are ranked between these. 
The major factors influencing the energy balance 


differences during the growing season are albedo, canopy 
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structure, snow accumulation, and growth characteristics. 
Short period variation in the observed differences reflects 
phenological differences among the six surfaces. Variations 
in the plots’ soil and air temperatures are a consequence of 
the aforementioned. 

Results are considered in relation to possible 
effects of rural land-use change on Central Alberta micro- 
climate. It is concluded that a change from farmland to 
woodlots will result in a decrease in albedo, and soil 
and surface temperatures, and an increase in canopy level 
air temperatures. Soil moisture levels will decline with 
increasing tree density, while an increase will occur if 
shelterbelts, low density woodlots, or shrub woodlots are 
used. Numerous small woodlots will result in an increase 
in sensible heat advection from surrounding farmland, de- 
pending on the woodlots' moisture supply. 

If present land use trends continue, and woodlots 
are replaced by farmland, albedo will increase, soil and 
surface temperatures will increase, while canopy level 
air temperatures will decrease. As a result, local con- 
vective activity may be enhanced, while soil moisture 
levels will decline, depending on the use of shelterbelts 
and/or snow fences. 

Differences between these and published results 
with respect to sensible and latent heat fluxes under 


woodlots, are due in large part to low soil moisture 
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levels generally existing at Ellerslie. Suggestions 
for future research should emphasize closer examination 
of: soil moisture and snow accumulation under wooded 
plots, canopy level air temperature, wind, water con- 
sumption of weed infested fallow vs tilled fallow, and 
surface reactions under unlimited or controlled soil 


moisture. 
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CHAPTER I 
INTRODUCTION 


Climate is changing on a global scale, and this 
idea has been generally accepted for many years. However, 
climate on a micro-scale is also undergoing change, and 
this can be traced (at least in part) to the influence 
of man. Although the effects of industrial and urban 
areaS on micro-climate have been widely recognized, man 
also affects his environment when the rural landscape is 
altered. By transforming forests into agricultural or 
pasture lands, and by stripping the vegetation in advance 
of mineral exploration and extraction activities. man has 
changed the rural micro-climate. This indirectly affects 
the watershed as well, since water yield and timing are 
dependent on the watershed'’s vegetation cover. 

Recent micro-climatological investigations of 
different rural land uses have become quite extensive. 
some of the research has concentrated on global and net 
radiation and albedo (Monteith and Szeicz, 1961; Stanhill et 
al., 1966;° Federer, 1968, 1971; Impens et al., 1970; 
Glover, 1972). Others have examined the flux parameters 
of the energy balance, i.e. net radiation, sensible heat, 
soil heat and evaporation (or latent heat), and subse- 
quently, potential and actual evapotranspiration (Fritschen, 
1966; Rosenberg, 1969b; Davies, 1972; Dilley and Shepherd, 
1972; Priestly and Taylor, 1972; Blad and Rosenberg, 1974). 
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While the above studies have many implications for 
both agriculture and watershed management, they are essen- 
tially examining static unique conditions. The energy 
' balance of a beanfield in Southern Ontario cannot be 
compared to that of alfalfa in the Great Plains because 
of the unique situation existing at each of the two lo- 
cations. Southern Ontario is a moist region, while the 
Great Plains are drier, and experience strong sensible 
heat advection during summer. Their respective soils 
have different structures and different moisture capac- 
ities. The surface topography and exposure of the plants 
may vary. There are latitudinal differences in radiation 
intensity and daylength. Most important of all, partic- 
ular weather events (eg rain) will not occur on these 
surfaces simultaneously, so the stresses they experience 
will be of unequal magnitude. Even studies comparing 
different natural vegetation covers in the same general 
location (Federer, 1968) suffer from the problem of un- 
equal exposure to climatic elements, differences in soil 
conditions, and diversity of plant communities. Muller's 
(1971) study of solar radiation in different coniferous 
stands also exposed this problem of comparison between 
different regions. Data from other studies could not be 
directly comparable because investigations had been 
carried out over a wide range of solar climates and 


because of largely unsatisfactory measures of stand 
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biomass which directly produce attenuation of the solar 
beam downward through the tree crowns. Another factor 
was differences in soil moisture, and the inability to 
separate evapotranspiration from drainage during moist 
periods (Satterlund, 1972). 

The word “static" was previously used to describe 
the above mentioned studies. This is due to the fact 
that the land use is constant throughout the research. 
With the many land use changes that have taken place in 
the past, there must have been changes in the micro- 
climate of these locations. Jeffrey (1964) has illus- 
trated several micro-climate changes that are caused by 
removal of forest cover. Muller (1966) showed that the 
abandonment of farms in New York State, and the subse- 
quent plant succession and reforestation of these lands, 
had affected the water yield of nearby watersheds. Ina 
review of recent research in watershed management, Leaf 
(1975) has indicated several micro-climate changes that 
can occur due to certain cutting patterns and practices, 
and reforestation. 

So, it can be concluded that land use change will 
affect micro-climate to some degree. Consequently, the 
change in micro-climate has some influence on various 
human endeavours such as agriculture, forestry, and land 
use planning. But research into the effects of land use 


change on micro-climate have produced more questions than 
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answers. Jeffrey (1964) lamented that "so little is 
confidently known of vegetation-water-climate inter- 
actions that a comprehensive list of research avenues 
_could be almost infinite." He included, among a list 

of priorities, an urgent need for research into albedo 
and energy exchange, and how they are affected by dif- 
ferent species, stand density, and structure. Baumgartner 
(1965) stated that energy balance studies of different 
plant covers were important since they could be used to 
predict the effects of different plant cover treatments. 
More recently, Satterlund (1972) echoed the need for more 
research on the effects of different plant species on 
energy balance. 

Additional uncertainties exist when trying to 
appraise the potential effects of altering the ratio of 
heating of the surface and of evaporation (the Bowen 
Ratio), especially if human activities, such as de- 
forestation and irrigation, are involved (Study of Man's 
Impact on Climate~-SMIC, 1971). 

Obviously, the need exists for research into the 
energy balance of different types of vegetation covers 
experiencing the same outside stresses. Under natural 
conditions, it is hard to imagine finding a pure conif- 
erous stand, a pure deciduous stand, agricultural crops 
or other plant covers sharing the same location, soil, 


and exposure. If this situation was common, comparison 
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of the effects of different vegetation covers on the 
energy balance could be accomplished easily, since all 
outside factors could be considered constant, hence 
isolating the effects of the vegetation from those of 
weather, soil, and topography. 

A unique and very suitable research site was 
available at the U. of Alberta climatological station 
at Ellerslie. This site was utilized for the purpose 
of investigating and quantifying gross differences in 
energy balance between the various forest, agricultural, 
grass, and fallow covers that are maintained at the 
station. Since topographic conditions and prevailing 
weather were identical for each plot, the differences 
in sensible and latent heat (evaporative) fluxes are 
directly comparable to one another. Thus, the potential 
effects of rural land use change on micro-climate, and 
subsequently, soil temperature and moisture, are 


examined. 
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CHAPTER II 


LITERATURE REVIEW 


Zek Introduction 


The following chapter includes a brief review of 
the background theory to be employed in this study. The 
energy and water balance equations basic to much of cli- 
matological research have been in use ever since instru- 
mentation and techniques were available to measure or 
derive their components. Besides examining these and 
related methods, this chapter deals with previous attempts 
at deriving individual components. Finally, there is a 
review of previous studies on the effects of land use 
change on micro-climate. Such studies may or may not 


have employed the methodology implied above. 


2.2 Radiation Balance 


lei Ver A Basic Physical Laws 
Physics has provided climatology with several 
fundamental laws on electromagnetic radiation. Wave- 


length is given by 
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The wavelength of maximum emission is related to the 
absolute temperature of the radiating body by Wien's 


Displacement Law: 


A max = 2.898 x ee cm °K (2.2) T = temperature 
. od sy 

These two equations yield both the wavelength and speed 
of radiation being emitted from the sun and the earth. 
From eq. (2.2), it is found that the earth's maximum 
emissions are at much longer wavelengths than the sun's 
because of the earth's lower surface temperature. 

Thestivux intensity of the radiation (1), is#a 
function of the fourth power of its absolute temperature, 


as given by the Stefan-Bolzman Law: 
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Eq. (2.3) is important when considering the relationship 
between the longwave radiation balance and surface temper- 


ature. 


24252 Radiation Balance Equation 


The radiation balance at the earth's surface is 


expressed as follows: 
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Rn = (Qtq) (l-o ) + In (2.4) where 
Rn = net balance of 
all-wave radia- 
tion) (Wm *) 
Qtq = global solar 
radiation (Wm) 
a = surface albedo 


In = net longwave 


radiation (Wm7*) 


This applies if both the atmosphere and the earth's surface 
are considered. 

Global solar radiation is made up of two com- 
ponents, direct-beam solar radiation (Q), and sky diffuse 
radiation (q)./ Of the solar radiation incident on the 
top of the atmosphere, an average of 22% reaches the 
earth's surface as diffuse, and 31% as direct-beam solar 
radiation (Sellers, 1965). Cloud cover affects the ratio 
of direct to diffuse by proportionately increasing the 
diffuse component and reducing the direct component. The 
former is almost independent of solar elevation. However, 
Since complete diffusion of light can be attained in nature 
only with a low overcast cloud cover and very low sun alti- 
tudes, measurements on most cloudy days will include some 
direct-beam radiation. 

Beer's Law describes the reduction of flux density 


as an exponential function of the depth of the homogeneous 


{* 


- 


anys? 2? ae 260 gilt: Dxts & 


c. : | . 
‘""}. noiveitex “Ss 


gvedian sortie = 2 


“pat) sink Selhins. 


Ao ata 


nelap insets “a ge 


Syewqpel tarr = ra | . 


retes nares ert (Seaphg ty Lhgee: i 
no. oa <. F. : . " SeCAtenss” ave ; 

=! F : =e 

gina Teiohy “9% 


~ =o Oat7, to. gi SHRM Bi. poltsitss 

: j 
wei PiSb (gan Sz || (oa ghebatbabe in blag. dah Some h.: _ Sheet 
ee 


7 " r oe" J ae “9 ian av 
GW: Se Jhebiocs no tle te . Se) 


ha iT : ' : 
é aie. doetidos YES Lo avi TEs 3 voligdanit: ‘oat ee 


aeice peed Ssh... as Sie bis »eenttth ss? Stabe a’ rae 40 
| (2304 veqetle®) nelvaibat 


ahtact: eft aifgatie sé¥00 Biol o 
ee wikeasitodi vista: Proqesy ye seuT Tih or toot 
FATE en}. aprkouh>: sofas: teramogMos 9¢ 


ae 4 
can 
lt ——— 
- mOaeEvEls les - 16 Prebheqabul taonle el’ 


tab iin ei ned 3 Zitat So neigvizib- aan 
‘ ‘a 
sali 


absorbing medium: 
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The extinction coefficient depends on the quantities and 
characteristics of the major absorbing and scattering 
materials within the medium, such as gases, water drop- 
lets, dust, or other suspended material. The situation 

is further complicated if the medium is a plant canopy, 
Since its architecture is rarely completely homogeneous. 
Besides this, the variation of species and their seasonal 
changes (such as the shedding of leaves by deciduous hard- 
woods) affect the canopy's geometry, and thus, the radiation 
balance of that surface. 

Changing weather conditions must also be considered 
when examining the radiation input available at different 
locations within a plant canopy. On cloudy days, the 
higher diffuse component of incoming radiation is instru- 
mental in creating a higher ratio of below canopy/above 


canopy radiation (Cohen, 1975). The absolute value of 
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radiative energy reaching a forest floor may not vary 
Significantly from cloudy to clear days (since the higher 
ratios occur during periods of lower energy input), there- 
‘by creating a daily radiation input that is much less 
variable than that of an open field. An open field 

has no overhead canopy cover that can intercept much of 
the direct beam component. Obviously, this must affect 
the longwave and energy balances as well. 

Surface albedo (a ) represents the fraction of 
incoming radiation reflected by the earth's surface upward 
into the atmosphere. Monteith (1959) showed that most agri- 
cultural crops which completely shade the ground have a mean 
albedo of 0.26, while most grass surfaces have an only 
slightly lower albedo of 0.25. Both values aronen remark- 
able consistency throughout the growing season. There will 
be some variation over the course of a clear day with a high 
direct-beam component, due to the roughness of the ground 
and canopy surfaces, the changes in sun angles, and the 
presence or absence of moisture. 

Albedo values for a variety of natural surfaces 
are shown in Table 2.1. In general, vegetative surfaces 
have higher albedos than bare ground. Albedo of forests 
(especially needleleaf forest) is, in most cases, slightly 
lower than both surfaces. Budyko (1974, p. 13) explains 
that this phenomenon is probably a result of better con- 


ditions of absorption within the limits of deep vegetation 
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cover, where there is an increased probability that solar 
radiation fluxes having penetrated the vegetation will be 
absorbed after their first reflection from elements of the 
vegetation. However, forests show marked albedo differences 
between species, and this could be a result of stand height 
and density, or reflectivity of individual leaves and 
needles. 

Although moisture content and solar elevation are 
major factors in albedo variance, prevailing weather is 
not, as the difference between mean albedo values of clear 
days and all days is negligible (Polavarapu, 1970). Albedo 
measurements taken just above the canopy give significantly 
higher values than airborne measurements, because of atmos- 
pheric attenuation of the reflection component, hence, the 
above canopy measurements are considered to be more pert- 
inent on a microscale (Rouse, 1965). 

Since a is the reflective component of global solar 
radiation, then (1 -, ) is considered the absorbable com- 
ponent, which makes up the largest proportion of the net 
radiation balance. This is contained in the short wave- 
lengths from 0.15 - 4.0 microns. 

The net. longwave component, unlike the shortwave 


radiation, is readily absorbed in the atmosphere at all 


iederailed review of albedo for various coniferous 
species can be found in a review by Jarvis et al. (1976). 
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but a narrow band of wavelengths, 8.5 - 11.0 microns, 
usually referred to as the atmospheric window, Longwave 
flux emissions can be related to temperature by eq. (2.3). 
Thessun issassumed to bewaablack body (e =, 1.0),  thateis, 
one which absorbs and emits all the radiation falling 
upon its surface. Most of the earth's natural surfaces 
are considered grey bodies with emissivities somewhat 
less than 1.0. Most forests and grasses, for example, 
have an emissivity of 0.90 (individual leaves may be 
higher), while values for most bare soil surfaces range 
from 0.89 - 0.98 depending on moisture content (Sellers, 
1965). 

Longwave energy emitted and absorbed by the 
atmosphere without clouds present is, in general, deter- 
mined by the presence of water vapour and carbon dioxide, 
and much less on ozone. The atmosphere actually emits 
longwave radiation both to the earth's surface and to 
space, but the downward flux usually originates ina 
lower and warmer atmosphere than does the upward flux, so 
the atmosphere will emit more energy to the earth's sur- 
face than to space. Therefore the longwave balance will 
usually be negative, and this, combined with the short- 
wave balance (which is positive during the day and zero 
at night) creates a net all-wave balance which is positive 
during the day and negative at night, depending on the sun'‘s 


elevation and the degree of cloudiness. From this, it is 
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evident that radiation income at the surface will also be 
affected by location, topography, and time of year, de- 
spite the fact that the sun emits radiation at a constant 


2nin + @lnOShLy. mins). 


rate of approximately 1,359.74 Wm 
Lambert's Cosine Law defines the radiation emitted as 
varying with the cosine of the angle between the normal 

to a surface and the direction of radiation. Therefore, 
an inclined surface oriented towards the sun will receive 
greater radiation intensity from the sun, and this has a 
large effect on the micro-climate of that surface. 

As a values for most agricultural crops are 
relatively constant (when there is full canopy coverage), 
so are the values of Rn. Freely transpiring vegetation 
does not have large differences in surface temperature. 
Therefore, if crops are in similar environmental condi- 
tions, small differences in longwave emission would result 
(Monteith, 1959). However, a change in albedo can be ex- 
pected to change the surface temperature, its longwave 


emission, and its net radiation. 


Pepe: Recent Work--Vegetation Transmission, Regression 
Analysis, and Variations in Albedo 


Several recent studies have focused on the effects 
of vegetation on the radiation balance. There were differ- 
ent objectives available for workers to pursue. Little 


work had previously been done on radiation distribution 
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in forests, and there were several unanswered questions. 
The spatial variation of direct, diffuse, and net radiation 
within plant canopies was unknown, and researchers desired 
to find this out in order to determine how much instru- 
mentation would be necessary to measure the eae avail- 
able within the total plant area for Snowmelt, evaporation, 
evapotranspiration, and photosynthesis. Other questions 
were: What was causing the spatial variation of radiation 
within the same plant canopy? Could this spatial variation 
be approximated statistically? How much of this radiation 
was being reflected? Was there a temporal] variation? 
Could water stress be a factor in net radiation or albedo? 
Also, Since global radiation measurements were more readily 
available than net radiation data, could global be used to 
predict net? 

Within the last decade, there have been attempts 
to resolve these questions. Anderson (1970) examined the 
spatial variation of solar radiation within forest and 
crop canopies at Madingly Wood, Zngland. She found that 
while the fraction of diffuse radiation was not much 
affected by weather and showed little spatial variation, 
the fraction of direct beam radiation transmitted varied 
as solar altitude changed. There was also considerable 
spatial and day to day variation because of local varia- 
tion of canopy structure and cloud duration. Reifsnyder 


et al. (1971) found that this extreme spatial variation 
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necessitated the use of 10 instruments to determine a 

full day average value under a pine canopy at New Haven, 
Conn., and many times that number would be required to 
estimate instantaneous radiation. For hardwoods, though, 
only one is needed for a full day average. Apparently, 

the different canopy structures of hardwoods and conifers 
yield different direct beam radiation profiles within them. 
While an exponential-extinction law similar to Beer's Law 
appeared to fit the coniferous canopy, a constant-ratio 
law sufficed for hardwoods. Therefore, while the conif- 
erous cover was acting as a uniformly scattering medium, 
the hardwood canopy acted more like a layer of horizontally 
flat leaves. Since the extinction coefficient of Beer's 
Law is a constant characterizing the optical opacity of 
the canopy, the difference in canopy structure is all 
important. This became a factor in the difference in 
instrument requirements for each canopy. 

Federer (1968) had found a range of 0.03 in the 
spatial variation of albedos in hardwoods, while for several 
vegetation covers, including conifers and herbaceous plants, 
the range was 0.13. The surface temperature range for 
hardwoods was 1 - 5°C, while a larger range, 5 - 11°C, 
existed for the other covers. These two factors were 
seen as the main cause of spatial variation in net radiation, 


Z 


which was 20.92 Wm tigre G0.03 sy, alee for hardwoods, and 


min” + (0, aL aliy, min7*) for the other covers except 


76.70 Wm 
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in winter. However, a more controversial conclusion 
arrived at was that one or two measurements of net radia- 
tion within each cover type were sufficient for calculating 
mean net radiation values, if the instrument was placed 
high enough to see the crowns of a number of trees. 
Reifsnyder et al- had disagreed, and Impens et al. (1970) 
later found that during sunny conditions with a high 
direct beam component, extreme values were generally 
much more frequent than mean values in several crop can- 
opies. This suggested that the determination of mean 
net radiation by few measurements may be grossly inade- 
quate. They concluded that several dozen instruments 
might be needed to determine even 12 hour totals of net 
radiation with reasonable accuracy. 

Besides this spatial variation, temporal variations 
of a diurnal and seasonal nature, were found for albedo. 
This variation is important because of its implications 
for net radiation. Kalma and Badham (1972) found that 
seasonal changes in the albedo of grasses were very likely 
due to changes in the surface characteristics of the vege- 
tation. On a diurnal scale, however, a general increase 
was observed with increasing zenith angle, and this must 
be due to the fact that natural surfaces are not perfect 
diffusers, and specular reflection becomes more important 
at lower sun elevations. This situation would only occur 


under periods with a high direct beam component as opposed 
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to diffuse. Under overcast conditions, the evenness of 
the light would not allow the different vegetation prop- 
erties to assert their respective influences on albedo. 
Another factor that appeared to affect albedo was 
water stress. Kalma and Badham (1972) found that wetting 
of the underlying soil caused a decrease in albedo, while 
under periods of water stress, an increase in albedo oc- 
curred, leading to a decrease in surface temperature 
(Kalma, 1972). Similar results were obtained by Glover 
(1972), although his findings were partially based on 
visual observation. However, he used regression analysis 
of net vs. global radiation to determine the albedo change 
between stressed and unstressed vegetation, and found an 
increase of 0.12 during the water stress period. 
Regression analysis was, by this time, widely 
used to estimate net radiation from global, because of 
its simplicity, and the resulting high correlations. This 


linear relationship can be expressed as follows: 


Hf 


Rn = a(Qtq)+bd C270) a 
bd 


slope 


{I 


intercept 


The slope and intercept terms are components of the equation 
of the line. Eq. (2.6) is generally applied to grass or 
crop surfaces, to day-time net radiation, and to the snow- 
free season. Accuracy is acceptable, although there may 

be problems with the intercept due to lack of data near 


the origin (sunrise and sunset periods when global is 
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near zero, and Rn is often negative). This error can be 
reduced, however, by applying daily or hourly radiation 
values. 

Attempts had previously been made to improve 
eq. (2.6), the most noteworthy being the heating co- 
efficient (Monteith and Szeicz, 1962). Idso (1968) 
claimed that this was invalid because it considered net 
longwave radiation loss to be a linear function of net 
radiation, which is a physically unrealistic assumption. 
Gay (1971) went further, when he stated that (Qtq) was 
being autocorrelated, and that regression analysis was 
producing useless results. However, other studies such 
as Glover's, continued to use eq. (2.6) because it was 
fulfilling its purpose. Albedo could be inserted into 
the regression model to produce different regression 
slopes, hence, different results for net radiation, be- 
cause of the variation in global radiation absorption of 


different surfaces. 


eis) Energy Balance 


ra eile Energy Balance Equation 
The energy balance at the earth's surface is 


expressed as follows: 


Rn = H+LE+Gst+ Y+M C273) H = sensible heat flux 
(wm7*) 
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a 


Le= ‘latentrhneat of 


vapourization 


(cal. em) 


ty 
" 


quantity of water 
evaporated (gm cm 
Gs= soil heat flux (Wm 
Y = photosynthesis (Wm 
M = miscellaneous ex- 


changes (Wm7*) 


The L and E terms combine to make up the flux of latent 
heat (units are cal om” which are converted to wm“), 
M includes energy exchanges which are due to metabolic 
activity and the storage of heat in the plant tissue or 
the volume of the canopy. 

Eq. (2.7) is an approximate term and does not 
include effects of snowmelt, which would affect the right 


é eye heat trans- 


side of the equation by about 6,970 Wm 
fer by precipitation, which would affect Gs, wind effects, 
which are an energy source and would be included on the 

left side of the equation, ane other biological influences 
such as oxidation of substances (e.g. by fire). Also, since 
the sum of Y and M terms is usually smaller than the experi- 


mental error in measurement of the major components, these 


are generally ignored so that 


Rn = H+LEt+Gs (eae) 
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adequately describes the energy budget of the earth's 
surface (Rosenberg, 1974). 

At night, with no shortwave radiation available, 
Rn equals Ln, which is usually negative since the earth 
is radiating more heat than it receives from cloud trans- 


missions. Therefore, 
In = H+LE+Gs Wak) 


The four components of eq. (2.8) can be directly 
measured or calculated from other parameters, such as 
temperature and vapour pressure (Sellers, 1965). Rn and 
Gs can be directly measured. H is calculated using the 


following equation: 


H = 60 paCp ah aT/aZ (2.10) pa = density of air 
Gives ax 10a. gm om7) 
Cp = specific heat of 
air at constant 
pressure (0.24 cal 
em? deg+) 
oh = molecular diffusiv- 
bir) Rene tsw Geo Swope? J eleheiey 
(0.16-0.24 om* sec”“) 
aT/aZ = temperature gradient 


within layer of 2 


thickness 
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This applies from the surface to the top of the laminar 
sublayer. Above that, the airflow becomes turbulent and 
oh is replaced by Kh, the coefficient of turbulent ex- 


2a 


change, which may vary from near zero to 10,000 cm 
& can then be calculated as the residual term of the energy 
balance equation. On the other hand, E can be measured 
using lysimeters, atmometers (porous ceramic or porous 
paper evaporating surfaces connected to a water reservoir), 
or evaporation pans. In that case, results for H can be 
determined as the residual term of the equation. 
"Evapotranspiration"(ET) is used to describe the 
total process of water transfer into the atmosphere from 
water and vegetated land surfaces. Soil evaporation (as 
well as evaporation of intercepted rainwater) and plant 
transpiration occur simultaneously in nature, and since 
the two terms cannot be distinguished easily, their com- 
bined effects are lumped into this term. ET can be esti- 


mated hydrologically, using a form of water balance 


equation: 


rainfall + irrigation - runoff - A PW -percolation 


= BT (Omit) 


where APW is the change in volume of water stored in the 
soil during some specified time period. There are problems 
with this because of measurement difficulties, as well as 


theoretical difficulties with percolation andA PW. However, 
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due to the climate and soil conditions of Central Alberta, 
runoff and percolation can be neglected. The local runoff 
hazard is low, due to the low precipitation and generally 
flat topography. Also, during the growing season, the 

soils of this area normally contain much less moisture 

than their moisture holding capacity. Therefore, the chances 
of deep percolation are also negligible (Verma, 1968). Since 


there is no irrigation here, eq. (2.11) becomes 


rainfall - APW = ET Ziel) 


Pare tes The Case for the Water Budget Approach 


In order to complete the calculations for the 
energy balance, either H or LE had to be determined by 
whatever means available. Eq. (2.10) would be very 
difficult to use because of the Kh term, which is ex- 
tremely variable. The temperature gradient, while easy 
to measure, would not be valid in this situation because 
OL sunemcharacteristics of the field) site to be used in 
this study. The plots range in area from 85 - 500 m”, 
and with such a small fetch, advective effects will be 
prominent for any temperature measurement taken above the 


plot surface. Therefore, the LE term will be calculated 


while H becomes the residual in eq. (2.7). 
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It now becomes necessary to decide the manner in 
which LE can be approached. Only the water balance (or 
water budget) method has been described so far, but there 
are many others, including the Bowen Ratio, the resistance 
approach, and the Penman and van Bavel combination methods. 

The Penman method involves the use of several 
empirical constants, but also requires knowledge of vapour 
pressures, sunshine duration, a , Rn, wind speed, and mean 
air temperature. The van Bavel method depends on the esti- 
mation of surface roughness. The resistance approach re- 
quires knowledge of air and leaf stomatal resistances, 
since these have profound influence on sensible heat trans- 
port and leaf transpiration respectively. This is an anal- 
ogy of Ohm's Law, since the latent flux is being regarded 
as a potential difference (in this case, T and vapour 
pressure gradients) divided by a resistance (Davies, 1972). 

Other methods, including the aerodynamic approach 
and the eddy correlation technique, have been hindered by 
instrumentation difficulties. An example of this is McBean's 
(1968) study of turbulent fluxes within a forest. Results of 
some recent studies appear to favor the Bowen Ratio, because 
the resistance approach has not been perfected yet, the 
Penman combination method has been shown to seriously under- 
estimate water loss under conditions of strong sensible heat 
advection, and the Van Bavel method was very sensitive to 


windiness (Rosenberg, 1969a, 1969b, Rose et al., 1972). It 
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thus becomes appropriate at this time to briefly examine 
the Bowen Ratio. 

Given the energy balance equation (eq. 2.7), and 
assuming soil moisture is not a limiting factor, so that 
Gs can be neglected since it constitutes only a small 
fraction of Rn, H and LE can be expressed as functions 
of the exchange coefficients Kh (for sensible heat) and 


Kv (for latent heat) respectively: 


Hite SC pie, iKhieeCan/aa)) tae 203) R 


pressure (mb) 
LE = aba ee Kv (de/dz ) (2.14) er = ratio of mole 
weights of moist 


airsto dry air 


i 


de/dZ = vapour pressure 
gradient within 
layer of Z thick- 
ness 


oa = density of air 


Other symbols have been previously defined. Subsequently, 


the Bowen Ratio (8) can be calculated thus: 


Wp ee Oe 
B = TR = a (aT /de) (2.15) 
Also; meicons Rn Ge beim Ges i knives 7 


An important assumption in simplifying the equation is 


Kh = Kv. 
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8 can vary between +o for a dry surface (LE = 0) 
to -l for an evaporating wet surface for which all the 
energy necessary for evaporation is supplied by convection. 
In other words, there is no net radiation gain so that 
H = -LE, and energy is directed towards the surface. 8 
may be very close to zero when LE is at a maximum, and 
LE = Rn-G. In the case of sensible heat advection, g is 
preater than or equal to 1, and LE is greater than Rn-Gs. 
Fritschen (1966) used the Bowen Ratio method on 
Sixsoltrerenteirrieaven crops ald tound that El erates 
were 1,.0-1.8 times that of Rn, indicating that large 
amounts of energy were being extracted from the air mass. 
He found the Bowen Ratio to be good for short period obser- 
vations, and one of the few methods that would yield valid 


results when crops were subjected to heavy and frequent 


irrigation, or where high water tables were present. 
Denmead and McIlroy (1970) compared evaporation calculated 
from the Bowen Ratio to lysimeter measurements and found 
the results to be accurate within +0.1 mm Tiana They also 
confirmed that the assumption of equality between the eddy 
diffusivities of heat and water vapour, upon which the Ratio 
calculations could be simplified, was valid (at least over 
the range of instabilities examined in the study). 

Fuchs and Tanner (1970) found a way to reduce the 


effect of a potential measurement error on the Ratio. 


Apparently, wet and dry-bulb temperature measurement 
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errors uweres ToundmaLo Deudirect ly propor.tionalatouthesRatio. 
The smaller the temperature gradient between the two measure- 
ment levels, the greater the potential error, which would be- 
come infinity at isothermy. If the temperature of the two 
air streams would be equalized, using differential psychro- 
meters, the measurement error would become independent of 
the calculated Ratio. 

Rouse and Wilson (1972) actually used the Bowen 
Ratio as the standard with which to compare another method's 
accuracy. This was somewhat premature, although it probably 
did not affect the conclusions they drew on the usefulness 
of the method examined, that being the water budget approach. 
It was found that this was reasonably accurate for long term 
periods, but could not give good daily estimates of ET, due 
to unavoidable large potential errors in measuring soil 
moisture. The use of neutron scattering devices was found 
to be more accurate than the gravimetric method, but ona 
day to day basis, this could not be relied upon for soil 
moisture measurements. 

Meanwhile, Dilley and Shepherd (1972) tried using 
a formula of the combination type, and the class A evapor- 
ation pan. These were compared with lysimeters, and both 
of them correlated well for estimating evaporation over 
irrigated grass and soil. Over a more variable and water 
sensitive surface, such as potatoes, correlations were some- 


what lower, but still satisfactory. However, there was high 
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local advection, and considering Rosenberg's findings 
(1969a, 1969b), and the local situation with respect to 

the siting of lysimeters, these correlations must be con- 
Sidered suspect. Davies (1972) also used the combination 
model for actual evaporation calculations over a non- 
irrigated beanfield. He found that the equilibrium evapor- 
ation results were closer to actual evaporation than to 
potential evaporation. = Since advection was not occurring 
at the site, the results can be considered valid. 

In a study examining non-steady-state conditions, 
when advection is prominent, Rose et al. (1972) found that 
the rate of energy stored within a crop canopy (V) was im- 
Pontans Isunisvsituatione it *was*irrefular “in’character; 
and often as large as Gs. Therefore, it may not be justifiable 
to neglect this in the energy balance equation. This could 
affect not only the combination model, but the Bowen Ratio 
as well. 

Blad and Rosenberg (1974), like Rouse and Wilson, 
had used the Bowen Ratio as a field standard, but it was 
first calibrated against a lysimeter, and underestimation 
of ET by about 20% was found during advective periods. This 


is in agreement with Rose's and Rosenberg's findings. It 


equilibrium and potential evaporation are special cases, 
requiring modification of the combination model. The equi- 
librium case occurs if the evaporating surface is located 
downwind from a long fetch over a moist surface. However, 
if the air next to the surface is saturated, this becomes 
the potential case (Davies, 1972). 
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was suggested tnat under these conditions, the heat and 
water vapour eddy diffusivities were not equal, and as 
previously mentioned, their equality is a factor in simpli- 
fying the calculations of the Ratio. Adjusting the ratio 
of these eddies to a value favouring the heat eddy would be 
sufficient to correct the observed underestimation. Ob- 
viously, these eddies would require further study, and 
until that time, the Bowen Ratio could not be considered 

a "universal" standard with which to compare other methods. 
That role is still occupied by the lysimeter. 

It becomes apparent that the Bowen Ratio will suffer 
from the same problem as the other approaches--advective 
effects caused by the small fetch of the plots at the field 
site. Due to the unavailability of an instrument as large 
and expensive as a lysimeter (and at least seven would have 
been necessary), the water budget becomes the most feasible 
method available. It is not significantly affected by the 
size of the plots, and equipment is available for the fre- 
quent soil moisture measurements that the method requires. 
The H/LE (or g) values, however, can be used for analysis 
purposes, Since they allow for quick comparisons of energy 
utilization to be made between surfaces. g values will be 


calculated from the final results for H and LE. 
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2.4 Effects of Land Use Change on Micro- and Macro-Climate 


OF a) Albedo and Energy Balance 


Stemming from previous discussion of items of the 
energy budget and the factors influencing it, one can now 
appreciate the significance of the effect of a land use 
change on the various items. A brief return to Table 2.1 
reveals that deforestation of deciduous trees in favour of 
agriculture or pasture may increase albedo as much as 10- 
14% because of greater bare soil exposure, while larger 
increases would occur if the original cover was bone feres. = 
Similar increases occur for revegetation of bare ground with 
grass. Deciduous reforestation of a bare area (as well as 
the reverse process) would have little short term effect 
on albedo. However, if the soil was to remain bare for 
long, and become a desert with a very dry surface crust, 
its albedo would sharply increase (SMIC, 1971; Satterlund, 
1972) 3 

A higher albedo will reduce the absorbable ebine 
ponent of eq. (2.4), thereby reducing the amount of energy 
available for heating the surface or for evapotranspiration. 
If the reverse land use change occurred, the resulting lower 


albedo would theoretically yield a higher surface temperature, 


lone example is Denmead's (1969) study of wheat and pine. 
He found that wheat's higher albedo led to lower Rn values 
because of greater exposure of bare ground in the wheat 
field. Unfortunately, the two sites were not examined 
Simultaneously. 
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depending on the soil moisture supply and other soil 
characteristics. It is important to note that the ground's 
surface layer acts as a heat reservoir--absorbing part of 
the excess heat of midday and summer, and returning it 
during the night and winter to the atmosphere. This stor- 
age lessens temperature extremes and acts as a climatic 
stabilizer, but of more short term importance is its abil- 
ity to radiate heat quickly. This creates steep atmospheric 
lapse rates, which can easily lead to convection and heavy 
tbumpopulence, 9.50 iisicenotesurprisine to find that higher 
values for H and lower LE values are common for bare soil 
(Satterlund, 1972). Meanwhile, vegetated surfaces with 
their lower radiation absorption rates (due to higher al- 
bedo) and better insulating qualities, show lower heat 
emissions and greater storage. This creates a stable lapse 
rate, and at night, inversions are more likely to occur 
over this surface, than over bare soil. Baumgartner's 
(1965) findings would seem to support this view. He com- 
pared the energy balances of several different surfaces 
in Germany, and found that bare soil had the highest soil 
heat flux rate. (Table 2.2). However, bare soil has the 
smallest Gs/Rn ratio, which could indicate that a land use 
change will not result in a linear change in either soil 
heat flux, or the other parameters of the energy balance. 
Another aspect of land use change is plant species 


and stand density. Dilley and Shepherd (1972) cite seasonal 
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change in crop structure, change in species combination, 
as well as plant physiology, surface geometry, and local 
environment as factors in land use effects on climate. 
Fritschen (1966) claimed that higher stand density leads 
to higher LE/Rn ratios because the reduced bare soil ex- 
posure provides greater insulation. This view was con- 
firmed by Satterlund (1972). However, Timmons et al (1966) 
found that corn evapotranspiration rates were not signifi- 
cantly affected by stand density, though it was noted that 
higher stands usually (but not always) had the highest 
evapotranspiration rates. The water budget approach was 
used in the Timmons study, with a sufficiently long time 
span employed between soil moisture measurements. 

In a review of recent literature on coniferous 
forests, Jarvis et al (1976) confirmed Dilley and Shepherd's 
views, citing cloud cover and canopy wetness as two important 
local environment factors. For most sites, irrespective of 
species, the daytime 8 of a dry canopy varied from 0.1 to 
1.5, and when the canopy was wet with rain or dew, between 
-0.7 and 0.4. However, certain locations show much larger 
values, exceeding those for field crops, and in some cases 
by a very large margin. It has been found that 8 values 
depend particularly on the stomatal resistance, but also 
on the general climatic and local weather conditions at 
the site. For continental climates such as Ellerslie's, 


higher 8&values are expected, as opposed to oceanic sites. 
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Jarvis et al. also quote values for LE/Rn at Fetteresso, U.K. 
When the canopy was wet, values were close to unity, since 
most of the water was used for evaporation, but under dry 
conditions, values ranged from 0.25 to 0.5. The authors 
conclude that the evaporation from dry foliage is only about 
5% of that from wet foliage exposed to the same weather be- 
cause the diffusion resistance of the canopy is much greater 
than that for air under dry conditions. 

Rauner's (1976) review of studies involving deciduous 
forests reveals that LE dominates the picture, especially 
during the second part of summer. During this period, small 
forests less than a kilometer in diameter create an "oasis" 
effect which reduces H to zero or even to stable negative 
quantities. In other words, the surrounding countryside 
supplies the forest with heat. This is also common in ir- 
rigated lands, and around springs and other moist areas 
under dry weather and in arid regions. Even in humid areas 
when Rn is low, advected sensible heat may contribute appre- 
ciably to evapotranspiration (Satterlund, 1972). For larger 
forests, this effect is rather weak. However, a similar 
effect can still occur if warm dry winds prevail over larger 
wet areas, even for a short time period. The chinook has 
been shown to affect a white spruce - fir forest at Marmot 
Creek, Alberta, in this manner (Storr et al., quoted by 


Wilson, 1974). 
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Seasonal change in vegetation structure affects 
values in deciduous forests. During early spring as leaves 
appear, g changes drastically from near unity to near zero, 
and maintains this level until defoliation in autumn. Rauner 
also states that deciduous forests seem able to use soil 
water very effectively. The structure of their heat bal- 
ance depends less on soil water than herbaceous vegetation 
or crops under comparable conditions. 

Recent grassland studies were reviewed by Ripley 
and Redmann (1976). At Matador, Saskatchewan, a 1971 
experiment revealed that 8 values approached unity during 
late summer, but this was more due to lower Rn and LE values, 
since H was fairly steady. LE/Rn was 0.7 in late July and 
0.5 in August, with Gs accounting for about 7% of Rn. Ap- 
parently, a large part of the decrease in LE was due to 
increasing transpiration and transpiring leaf area (mainly 
the result of water stress and high leaf temperatures). 

The remainder was attributed to the decrease in evaporation 
from the soil due to the drying of the upper layers. De- 
spite the fact that much of the total root biomass is above 
30 cm, the steadying of 8 near unity during drought con- 
ditions in late August indicates that the vegetation was 
able to use soil water from lower depths. Another explan- 
ation put forward by the authors is that stomatal closure 
and perhaps internal plant and soil resistances would limit 


evaporation during much of the day during the evaporative 
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demand period of late July and early August. These effects 
would be reduced in the cooler month of September and evap- 


oration might continue through the midday period. 


Bnergy balance changes can act as indicators of 
alterations of other climate parameters. Fritschen (1966) 
concluded that higher LE/Rn ratios indicated warmer air 
temperatures, while lower ratios indicated higher atmos- 
pheric vapour pressure. Higher ratios, especially those 
greater than 1.0, indicate that energy is being extracted 
from the air mass by the plants. On that basis, Fritschen's 
experimental results indicated to him that broad leaf plants 
use more water than grasses. He attributed this to their 
rougher and more rigid surface structure which would allow 
greater energy extraction from the air mass in the form of 
sensible heat. 

Table 2.3 illustrates a theoretical change in energy 
balance if forest is converted to agricultural use. This 


° * ° e e e -2 
assumes an incoming shortwave radiation intensity of 125 Wm -, 


and an outgoing infrared radiation level of 50 Wm", It is 
interesting to notice the difference between wet and dry 
arable surfaces. Apparently, the first priority of energy 
use is to evaporate water, hence a low g value. Once ET 

is satisfied, the remainder of the energy is used to heat 
the surface, and g increases. 


Natural grassland, which is found in semi-arid and 


prairie regions had the highest g value, but this is still 
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well below typical desert values of 2.0 - 6.0 or higher. 
Coniferous forests have higher sensible heat values than 
deciduous, due to its lower albedo. This will probably 
‘mean warmer temperatures during high radiation periods. 
During periods of low radiation, the surface's emissivity 
and insulating qualities determine temperature levels. 

Many scientists believe that the energy balances 
created by widespread deforestation could be enough to 
significantly alter circulation patterns, especially in 
tropical areas (SMIC, 1971). In the mid-latitudes, espe- 
cially the prairie regions which are only slightly wetter 
than deserts, Bowen Ratios become very important indicators 
Of=clumate trends. © fi 6 1s*preaver than 1.0. it indicates 
very high sensible heat levels, low availability of moisture, 
and a warm continental arid climate. 

At the present time, little replacement of temperate 
forests by arable lands is in progress. Indeed, there are 
many marginal farm areas that have seen a return to forest 
cover, such as the Piedmont area of Georgia (Brender, 1974), 
Central Manitoba (White, 1975), and Northern New York 
(Muller, 1966). However, urban expansion and mining ex- 
ploration may create a very different trend, with its own 


unique effects on climate. 
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2.4.2 The Hydrologic Cycle 


Land use change alters the surface cover, and so, 
has a profound effect on the various components of the 
hydrologic cycle. But many of the changes are not obvious, 
and in the final “analysis; it is difficult to conclude 
whether there will be an increase in water yield, or a de- 
crease, without first considering a wide range of factors. 

Deforestation is a fairly drastic land use change, 
whether conifers or deciduous trees were the original cover. 
So the question is asked--what changes will occur if an 
area's forest cover is removed? To answer that, each hydro- 
logic component, as well as the area's forest cover must be 
considered. Air temperature, wind, and absolute humidity 
are affected, too, and since these factors influence evapo- 
transpiration, they are discussed here. 

In general, deforestation results in increased low- 
level (low altitude) maximum temperatures of from 1 - 5°C 
in the warm season and from 1 - 2°C in the cold season. The 
warm season temperature change is dependent on forest type. 
In addition, minimum temperatures decrease 1 - 4°C in all 
seasons. Therefore, al - 9°C greater temperature range 
is produced on a monthly and annual basis (Changnon, 1972). 

Absolute humidity at and near the surface is re- 
duced 2 - 25%, depending on cover, forest type, and the 
new land use type. The amount of moisture added to the 


atmosphere by evapotranspiration is lessened, but since 
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low level winds in deforested areas are increased 2 - 4 
times, a greater rate of evaporation (and more wind erosion) 
from the soil occurs, partially offsetting the reduced evapo- 
transpiration. It has been suggested, that increase in 
drought frequency and a general drying of the climate in 
eastern Europe have been due largely to deforestation and 

to poor surface water management, which together led to a 
lowering of the soil moisture. Other controversial questions 
concerning land use and climate include deforestation's 
effect on rainfall production in the Mediterranean, and 
climate changes in semi-arid regions caused by cultivation 
of previously undisturbed grasslands (Changnon, 1972). 

These are macro-scale phenomena, and are beyond the scope 

of this study, but they are worth mentioning because they 
illustrate the potential magnitude of the problem. 

The components of the hydrologic cycle include 
interception,, throughfall, infiltration, prunoff, and evapo- 
transpiration. If deforestation occurs, interception losses 
will be reduced and throughfall will increase. If the area 
is on a slope, this would affect surface runoff and erosion. 
However, the new vegetation will be different in density 
and Senne arrangement so as to reduce the effects of wind 
while it is becoming established. These "replacement plants" 
may transmit or reflect more radiant energy from leaf sur- 
faces, or transfer absorbed heat to surroundings more ef- 


fectively, so that less is used in transpiration. Also, 
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they may have a shorter period of soil moisture use or a 
shallower or perhaps, less efficient root system, and thus 
extract less water from soil storage (Jeffrey, 1964). 

On the other hand, not all grasses have shallower 
root systems than trees or shrubs. Grass species vary 
from one to the other, just as do trees. Comparisons 
must be made--species vs species--for the kinds of sites 
under consideration. For example, the conifers have not 
been shown to use either more or less water than hardwood 
trees, with similar rooting depths (Jeffrey, 1964; Satterlund, 
1972). Trees and shrubs growing along stream banks and in 
areas of high water tables have been shown to use large 
amounts of water. Transpiration varies with age too, in- 
creasing up to a certain age, levelling off, and then de- 
creasing (Jeffrey, 1964). 

Infiltration will be greatly affected by age, stand 
density, and human activity. High infiltration rates occur 
in dense, old, and undisturbed stands. Highest rates occur 
in undisturbed forest while lowest rates occur in overgrazed 
pastures and in cropped fields. Rates will increase if the 
area is reforested, thereby reducing runoff. Consequently, 
streamflow leaving the drainage basin will be reduced 
(Muller, 1966). 

Time of year will be a factor when comparing the 
water use of different tree species in the same location. 


Satterlund (1972), stated that the seasonal pattern of 
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soil moisture removal is different for aspen and spruce in 
Colorado, or red pine and oak in Michigan, to name two 
examples. 

Despite the effects different plant species may 
cause, there is still one basic similarity among them-- 
their presence on the surface and their roots below. If 
the soil is bared, only a very shallow depth is dried since 
there is less energy available to transfer deep into the 
soil (depending on albedo and H values), as well as the 
lack ‘of ‘roots to cause water uptake (Satterlund, 1972: 
Muller, 1966). Farmers in prairie and semi-arid environ- 
ments often use fallow rotations (dry farming) as a water 
conservation practice because of this fact. 

Having examined each of the separate components, 
one finds that there are both positive and negative feed- 
backs from each when a land use change occurs. What is 
the net result of these feedbacks? Muller's (1966) study 
of abandoned New York State farmland which has returned to 
natural overgrown pasture, deciduous forest cover, or con- 
iferous plantations, reveals a decrease in water yield. 
Both in winter and summer, radiative energy and moisture 
are more readily available for evapotranspiration in these 
land uses. During winter, greater snow evaporation occurs 
in snow-covered crowns of conifers, because of favourable 
daytime radiation balance. In summer, deciduous and con- 


iferous trees absorb more radiation than open and overgrown 
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areas. Also, differential rates of evaporation and trans- 
piration during moisture deficits, and greater infiltration 
rates, will increase evapotranspiration loss and decrease 
water yield from the overgrown and forest covers. The 
greatest changes occurred when conifers were planted, as 
water yields decreased by up to 24.4 cm yr7t, Muller 
had to qualify his conclusions, however, since subsequent 
water balance and multiple regression tests were signifi- 
cantly affected by estimates of precipitation income and 
water yield outgo which may have been unrepresentative of 
watershed-wide conditions. Mean monthly temperatures from 
an open base station were used for evapotranspiration cal- 
culations, and with reforestation, rain gauge exposure 
changed. Therefore, it was admitted that land use change 
may not have been the exclusive factor. 

Taking these factors into account, this study will 
use mean weekly temperatures from all different cover sites 
when making evapotranspiration calculations. The open 
station data will be used for open area sites such as 
fallow and shrub, while data are available for the individual 
spruce and poplar plots. Only 1 rain gauge is used at the 
site, but since the plots are all within 500 metres of the 
rain gauge, significant variation is not expected. Also, 
as previously mentioned, runoff and percolation can be neg- 


lected under the local conditions studied. 
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INSTRUMENTATION AND EXPERIMENTAL PROCEDURE 
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The following presents a description of the field 
site and test plots to be examined during the course of 
this study. A detailed account of the type of instru- 
mentation used, and its deployment, is included. Cali- 
brations are listed in Appendix I. A detailed chronicle 
of the biological changes observed on each plot from 


Apri MoeOctober 1 1976,e1s found: in Chapter iV. 


B.2 The Site 


3.2.1 General Characteristics 


The site chosen for the investigation is the 
University of Alberta Farm at Ellerslie, Alberta (lat. 
S32 5 aN Lone. g Meee jolla Wy eabdutaens km. south. of 
Edmonton “secity core (Fig. 93.1), Located in the trans- 
ition zone between the Boreal and Prairie regions, the 
natural vegetation here is mostly trembling aspen (Popu- 
lus tremuloides Michx.) and white spruce (Picea glauca 
Voss). However, a large portion of the area is under 
pasture and mixed grain cultivation with barley being 


the most popular grain crop. 
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The area is underlain by shale, siltstone, and 
sandstone deposits from the upper Cretaceous period. 
Ground moraines, creek channels, and prairie mounds are 
evidence of previous glacial activity. Topography at 
Ellerslie is fairly level with slopes rarely exceeding 
2% (Fig. 3.2), except at river and creek channels. The 
experimental plots are on a NNW facing slope of about 1%. 

A cool continental climate exists here (Dfc on 
the Koppen system) with 2 mean summer temperature of 
13°C, and mean summer precipitation of 30 cm. The average 
growing season is 165 days, commencing April 25 and ending 
October 6. However, only 90 days are frost-free. June 
and July are the wettest months, but rainfall is generally 
low in intensity and well distributed throughout the 
growing season. There is a comparatively low variability 
in growing season precipitation here. Average variability 
is 20 - 25% from the mean, whereas elsewhere in Alberta, 
especially in the northwest and south, this variation 
exceeds 35% (Weir, 1971). There are about 110 - 120 
days with precipitation per year, mostly during summer. 
Mean annual soil moisture deficit is 10 cm and mean annual 
potential evapotranspiration is 51 cm, 5 cm greater than 
mean annual precipitation. Annual bright sunshine hours 
slightly exceed 2200. 

The annual march of the components of the energy 


balance for the continental mid-latitudes shows differences 
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in time of maximum flux intensity. A sample curve 
(Bara loo Reo ySnowlein rie. 4.38) illustrates a 
large time difference between LE and Rn maxima, which, 

in this case, is close to 6 weeks. The maximum outgo 

of heat for evaporation can be observed earlier or later 
than the maximum in Rn, depending on the annual march of 
precipitation and the associated change in soil moisture. 
If the LE maximum precedes the Rn maximum, then the H 
maximum is commonly observed after the Rn peak, and vice 
versa. For climates with wet summers, heavy cloudiness 
in summer reduces the Rn maximum (Fig. 3.3b). The LE 
max. occurs in August, while the winter period shows 

a Curve similar to Fig. 3.3a (Budyko, 1974). Although 
this curve is more typical of a mid-latitude summer mon- 
soon climate, there are similarities in annual precipi- 
tation distribution and latitude between Grossevichi and 
Ellerslie. Ellerslie's annual curve for any particular 
year could be similar to either of the above two figures, 
depending on summer precipitation levels. 

Soil climate is moderately cold cryoboreail. Mean 
annual soil temperature is 2 - 8°C, 120 - 220 days are 
above 5°C, there is a short thermal period of less than 
60 days with temperatures above 15°C, and there is a mild 
summer. Soil moisture regime is subhumid, with a 5 - 10 cm 


water deficit during the growing season (Fremlin, 1974). 
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FIG. 3.3: ANNUAL MARCH OF HEAT BALANCE 
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Soil types in the small plot area are Malmo clay 
loam and Ponoka loam. Malmo clay loam is developed from 
lacustrine parent material, generally underlain by till 
a depths varying from 51 - 102 cm below the surface. 
The upper horizons of the profile (Ap, Ah), are black 
(10 YR. 2/1 - 3/1) with a clay loam texture and coarse 
granular structure. Ponoka loam is coarser in texture, 
with a parent material of stratified silts and very fine 
sand clay loam. Till generally occurs at depths varying 
from 38 - 102 cm below the surface. The Ap and Ah hori- 
zons are of Similar colour, but the texture is loam rather 


than clay loam (Pawluk, 1964). 


3.2.2 Experimental Plots 


Six small plots and one forest site were employed 
during the studya(Pieee3.2).s sAudeciduoussplot (PS), 
330 m2 in area, consists of poplar trees approximately 
6>metres highmarlatcomemlasandae) Dom Outh of it As a 
coniferous plot (S) approximately the same area (Plate 
3.2). White spruce is the tree species occupying it. Trees 
are 4 - 6 m high and form a very dense canopy. A shrub 
DLO te (W)hpel5O mic consists of western snowberry less than 1m 
high, tall grasses, and some weeds, including Hawk's beard 
(genus Hieracium) and Canada thistle (Carduus arvensis) 
(Plate 3.0)55 ebanley “Bjeand ftablow plots <(F) ,9S5 m< 


each, are the only ones which suffer seasonal human 


es § mbaraomy Etienton 
Jeon aa eh voted mo Bot - 1 wont galley iets 
Hout ‘ons. (HA. aR) otddone ont to snosiied teqqe ea? 
desc Mes coudeey mot unto «le CINE + Ne #Y OL) 
\sierye? mt aeargoo ‘at anol efteet -ouutoutta telunsta 
cei view bua atic bebiitente to Lalvetent thgeeq 2 Atty 


aribyitev eAtqee t= 207 viiermes: ifit .maot yaks bras | 


=}for if bas is 347 Sree att woled mo SOL - Bt mort 
veitRt mot al syurnet scr dod “eetoe wilteia to ets enoz 
(OSE. hectieaty/ cea 


‘gaieabadiaame sae 

beyatqmy giaw shin saad sito few atoi¢ £Sinino ake 
(28) tole ana tnt ane A , (9.0 «ght ybure%ed? gataub 
efotenixorags e0s%t xeldeg Yo edetenns yapem at Sy OFE 
gat si) 2g wreoe, AGerk men ae of ejaft) git cotton 8 
siell) sare ontwoewt yfatamdxongge (8) tol Suotetinee 
egext #1 eniuquoro eéioeqc agit af? of oowaoe OfiAW «(S.E 
dintde fh saree sezceb wer fio fan at ao - 9 om 
sik nah toeer vatedvire exited xo eaitedos Se ORE <luh sole 
breed 2 owl aethyiga! yebeeW aque hae esseety Liat yitalll 
~~ pinistlosed bis (omesowedit wines) 
Sa 28 nae eae yotred «(C4€ erat) 
Holly geno. Mao ott ets ,fose 


53 


isible 


LS eva 


The screen 


Deciduous plot (PS), 


Onmtneomlette 


Plate 3.la. 


Piawe 3.1). 


PACE Biecs 


ineidesro splot. From lett ito rieht: ‘screen, 
soil moisture access tube (at author's feet), 
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Coniferous plot (S). The tower is visible 
at the center. 
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Plate 3%. 


Plate 3.4. 


Shrub plot (WwW). 


Barley (B) and fallow (F) roots (May, 1976). 


Barley stubble (background 
plot was seeded one week aft 


taken. 
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Poplar forest (P). 


Plate, 9.6. 
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Plate 3.7. 


Ellerslie small plots (July, 1976). 
to bottom: PS, S, W, B (with radiati 
ment) and F (with heavy weed cover). 
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attention. Stubble from the previous barley crop was 
removed in spring 1976 and the plot was re-seeded, 

creating a continuous cover slightly less than 1m _ high 
by June. The neighbouring plot (F) was kept fallow through- 
out the experiment, although periodic weeding was necessary 
since common weeds such as Ball Mustard (genus Brassica) 
and Stinkweed (Datura stramonium) invaded the plot (Plate 
3.4) A short grass cover (G) was maintained as part of 

a well-serviced meteorological installation south of the 
fallow plot (Plate 3.5) The actual measurement area was 
approximately 25m from F plot. Distances between the 
other small piots ranged from 5 - 13 m. 

In addition to these sites, a nearby poplar forest 
site (P) was used on a limited scale. The site is located 
approximately 300 m west of the small plots, and consists 
of poplars 14 — 18imetresthigh (Plate 3.6). This forest 
area is much older and larger than the small plots, and 
has a denser cover than the small poplar site. Table 3.1 
summarizes the above descriptions, and includes the plot 
symbols to be used throughout the text. Plate 3.7 shows 


an overview of five of the small plots. 


SRS} Instrumentation 
Joel) wRhadilablon 


Global solar radiation was measured with an Eppley 
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pyranometer (Plate 3.8), sensitive to shortwave radiation 
in the range 0.3 - 3.5,. The sensor is constructed with 
a circular multi-junction copper-constantan thermopile 
-with temperature compensation, housed in 2 optical glass 
domes. It has an approximate sensitivity of 9 microvolts 
per Wm <, Reflected solar radiation (albedo) was measured 
with a second Eppley pyranometer, held inverted approxi- 
mately 1 metre above the ground surface for G, B, F, and 

W plots, (Plate 3.9). For Pseand S plots, the pyranometer 
was attached to the top rung of an instrument tower (Plates 
3.10 - 11). The resulting voltage reading was recorded in 
an electrically powered Leeds and Northrup strip chart 
recorder (Plate 3.12a). Digital counters were used for 
both daily global and albedo measurements (Plate 3.12b). 
These provided daily voltage totals. The albedo recorder 
was housed in a white shelter located next to the plot 
being measured, while the global Honeywell recorder was 
housed in the lab building. 

Net radiation in the range of 0.3 - 60.01 was 
measured with two CSIRO net pyrradiometers. The sensing 
element is similar to that of the pyranometer, but since 
it is housed in plastic hemispheric domes, the longer waves 
are permitted to strike the sensor directly. The plastic 
hemispheres cover both sides of the element (Plate 3.9), 
so that incoming and outgoing radiation are measured sim- 


ultaneously. The domes were kept inflated by an air pump 


61 


bis: a (3 nok eget 
retsmonnnye sar fegeia & bee RE eam (Cie BLD Gator W 
agte ll) owot eer oder gor ent of batéatte saw 
ar ‘pwbrrose% arnt aed Batti apahiow dativest eit ME = OLE 
daotis. qhahe qustdem, Peue chest vevawaq-eliss fitoobe da 
40% wae aia afedfbi catiaia eatie obit) asbi10091 
GSS. og aitrsne maha Ciba ban Kelis vhbed Woe 
ees obad is ani ager syrah hee {lias bepivere scent 
jotc ott oF daeq Geteoch -vestade otldw 2 Wt ‘beauort 2aw 
2a sabre sete Eigeirgnoit addin ert ete shea SEES 
| sgntehivd def eit eb besvon 


saw 1 0,08 0X egret ett ot robtethet tet © as 


nieces of? Jpubtonmelbemyg fan ORTEO ow? dt iw Betoasem 
vohie tu ,vedemeeree odt'ro thdt of tal inte Re tromete 
nevaw wagnol eft ,seatoh cinemnelmed oitasla ni beaued wl St 
obteaiq of? Aftoowlh tomas edt oniats of bettfareg q ete 
AG.t walt) tnpaste od?! t0 table Atod xevoo eexeigelaen 
“<a bonsnnen: eum soksnther,yelostus Sou sniaianh tae oe 
eae einai | igptow Semob edT sy fewoonet ly 
5 


— a “2 - aa 


7 _ 
a 


62 


Eppley pyranometer. 


Elate 3787 


Plate 3.9 Albedo and net radiation measurements. 


Plate 3.10. Albedo measurement above deéciqduouseplot. 
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Instrument tower at S$ plot: 


Plate 3.11. 
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Platen3)1 2a" Albedo recorder. 
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Plate 3.12b. Digital counter. 
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that allowed a flow of dry air to be continuously main- 
tained. Silica gel was used as a drying agent (Plate 
3.13). The resulting net all-wave balance was recorded 
on a Rustrak battery-operated strip chart recorder (Plate 
3.14), housed in a white shelter similar to the afore- 


mentioned. 


Be. 2 2Oil) Heat hilux 


Soil heat flux was measured with a soil heat flux 
disk, similar to the one used by Fuchs and Tanner (1968). 
The disk is circular and quite small, and is inserted be- 
tween 1 and 3 cm below the soil surface (Plate 3.15). 
The upper and lower surfaces of the disk attain a temper- 
ature difference which is proportional to the soil heat 
flux. This temperature difference is sensed by a thermo- 
pile similar to those used for radiation measurements. 
The output was recorded on the same type of battery- 
powered recorder used for the net radiation output 


(Plate 3.14). 


3ro63 voll Moisture BulksDensity, Organics Matter 


Surface soil moisture was measured with a Nuclear 
Chicago surface mcisture probe. The probe has a Radium- 
Beryllium source which emits fast neutrons into a hemi- 


spherical volume of soil. When the neutrons are 
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Plate 3.14.,,Automatic recordin 


g of net radiation and 
soil heat flux. 


Plate 3.15. Soil heat flux disk (rotated). 
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sufficiently slowed by the hydrogen nuclei in soil 
moisture, the detectors will produce a pulse which is 
counted by the scalar (Plate 3.16). The number of counts 
per unit time (in this case, one minute) was divided by 
the standard count, which was accomplished using a granite 
block and paraffin standard (Plate 3.17). These were taken 
before and after each set of readings. The resulting ratios 
were converted to percent volume of soil moisture by cali- 
brating the probe against gravimetric measurements taken 
in April. Bulk density measurements of the soil allowed 
percent volume to be converted to cm water “ee soil, there- 
by facilitating data use in the water balance equation. 
Soil *moreturesateo. 2,905, 059 and 1.0 metres 
was measured with a Nuclear Chicago depth moisture probe 
(Plate 3.18). It has a source similar to the surface probe. 
Data were calculated by using the “ratio method" as above, 
with Alberta Research Council calibration standards being 
employed rather than gravimetric samples. One access 
tube for each small plot, and three for the forest plot 
were installed at the end of April, 1976. Measurements 
of surface and depth moisture levels were made each week 
from May 2 - October 7, 1976. 
Bulk density (Db) samples were collected with a 
Bull Core Sampler in June. The oven dry weight of the 
core samples was related to the core volumes to calculate 
dry density. Subsequently, porosity was calculated as 


follows: 
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Piatems:.16. Sin Seas moisture probe and scalar in use at 
plos. 


Plate 3.17. Standard count apparatus for surface probe. 
The probe is sitting in the paraffin which 
is placed on the granite block. 
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Plate 3;18; 


Depth moisture probe with scalar in 
at F plot. 
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Db 


POR% (1 - Sn eyes) x 100 Gar) 

The constant is the assumed particle density of 
soil. For surface horizons, this was lowered to 2.55g/cc. 
Results are listed jim Table 3%? 

Percent organic matter content was determined by 


ignition at 450°C. Results are listed in Table 3.3. 


3.3.4 Temperature and Relative Humidity 

Maximum and minimum temperatures were measured 
with Taylor or Casella max-min. thermometers. These were 
placed in Stevenson screens, which were installed approxi- 
mately 1.7 meters above the ground in PS, S, and P plots 
in May. Data for the other sites were provided by equip- 
ment in operation at the Department's meteorological in- 
stallation, which doubled as G site (Plate 3.19). 

Relative humidity at PS, 5, and P plots was measured 
with Negretti and Zambra or Casella thermohygrographs, lo- 
cated in the same screens as the thermometers. (Plates 
3.20, 3.1b). Wet and dry bulb thermometers were used 


for the open site measurements (Plate 3.19). 


Se Ste) Precipitation and Evaporation 
Precipitation was measured with a standard rain 


gauge located at the meteorological installation. Open 
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Blate 3.9. 


Maximum and minimum temperature and relative 


humidity measurements at G plot (the meteoro- 
logical installation). 
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Plate 


3.20. 


Thermohygrograph and maximum-minimum thermom- 
etereinescreen at Ps plot. 
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air evaporation was measured using a Class A evaporation 


pan (Plate 3.21). 


SS Deployment of Zquipment 


One net pyrradiometer and one soil heat flux disk 
were installed at G plot on May 22, 1976. These remained 
stationary throughout the measurement period, and so be- 
came the base station. The pyranometer (global) was also 
stationary. The other three radiation instruments were 
rotated among PSs sie, 5b, and@replote-sm During the May — 
July period, the net pyrradiometer and albedo instruments 
were moved twice weekly as a precaution against creating 
too large a time gap between readings for any individual 
plot. During August, this was changed to weekly. The 
soil heat flux disk was rotated weekly due to the sensi- 
tivity of the instrument. 

Surface soil moisture was measured within 24 hours 
of previous precipitation from May 22 - July 11, 1976. 
For the remainder of the program, surface readings fol- 
lowed the same schedule as the depth readings. 

Relative humidity was measured at 0900 hrs. MST. 
Spot measurements were effected at G plot, while points 
were taken from the thermohygrograph charts. These 
chart readings are probably not very accurate due to cali- 
bration difficulties. The data will be reviewed briefly 


in Chapter IV. 
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MS 


Evaporation pan (foreground) and rain gauge. 
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So: SIs Supplementary Measurements 


Additional data were obtained from a neighbouring 
large barley field, September 5 - 10, 1976 (Plate 3.22). 
~The barley had reached maturity, and had very few weeds. 
Net radiation, albedo, and soil heat flux were measured 
continuously (Plate 3.23). Surface soil moisture was 
measured September 5 and 10. Concurrent moisture data 
were obtained at B and G plots, and at a grass area mid- 


way between B plot and the large barley field. 
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Plate 3.22. Large barley field. P site is at the upper 
left. Instruments and recorder shelters are 
on the right. 
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Plate 3.23. Albedo and net radiation instruments on large 


barley field. PS and S plots are on the right. 


The lab building is on the left. 
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CHAPTER IV 


GENERAL CLIMATE AND FIELD OBSERVATIONS 


4.1 Introduction 


This chapter deals with general climate and field 
conditions obtained before and during the field program. 
Climate parameters to be discussed are air and soil 
temperatures, precipitation, evaporation, and sunshine 
hours. Current trends aré compared with previous data 
collected over the past 12 years at Ellerslie. Changes 
in soil moisture and vegetation characteristics during 
the growing season are described. 

Soil temperature data comparisons are for the four 
year period 1970 - 1973 instead of the 12 year period. 
This is due to the relatively incomplete stage of the 
tree plots’ canopy development before i970, and the lack 
of data for 1974. 

Relative humidity is briefly discussed in Section 


4.2.6. Data comparisons are for the study period only. 


4,2 Twelve Year Mean Climate Conditions 


4.2.1 Air Temperature 
Mean annual air temperature for Ellerslie was 
1.4°C. Monthly values ranged from -17.7°C in January to 


16°C in July. Mean monthly values below 0°C have been 


82 


ed ° 
7 ‘ 
; oe |) : 
> io ay uf > ae ave 
te in _ 
3 2 f ee ad eS ae 54 ; ‘ 
, a 


» j 
er ee ¥ 5 %, cs “ t 


prot has stamiie, Kenenog dtle elash. * 
margorg DLott ott gabaub bow eicred, 
Chem. Meee subse eum becusseth..of -9Fne: 
eninanve beus aobteroqart veolsaalannt 
 ptab_ avolverq a20 beso WTA. eet 
gepued? .ollavetie-te eamex St Ses RE, 
gcrdnpant Deena aee me ek. 
; bedfroend: or snd ed way wr 
agoy ait wok o1e enveltaqnoe ‘niab-oiwtarsame iiee 

olive tae SE agit, 10, aaweak, €Iet = ORCL hole 1887 
a to seats etelqnoont ufavitsrar ett of “esb at: ‘ghaT 
dnal wat dirs 1 ee ae 
si we ated To 

ngitsee al —— ef te hed, oh Yribknned ov tation | 
eine holzed youte aA t0% ous snogkaagmos stad 2.3.4 


=_ 


recorded from November to March (Fig. 4.1), and individual 
negative monthly values have occurred twice during April. 
Seasonal change appeared to be fairly regular, although 
deviations could be seen in March and September, both of 
which appeared slightly cooler than might have been antici- 
pated. 

Winter monthly means were highly variable. December 
had recorded values of -5.5°C and -21.3°C, while extremes 
for March were 0.7°C and -11.5°C. During the summer months, 
however, temperatures were more consistent. July, for ex- 
ample, had recorded extreme values of 14.0°C and 17.7°C. 

Since the available data cover a short time period, 
winter mean monthly temperatures should only serve as rough 
approximations, rather than useful mean values. Due to the 
wide range of individual values, a larger sample size is 
required before means can be calculated with confidence. 
Summer mean monthly temperatures are more reliable, how- 


ever, because their ranges of individual values are small. 


PotD Precipitation 


Definite wet and dry seasons are evident in the 
mean monthly precipitation curve (Fig. 4.1). Mean annual 
precipitation was 454.7 mm, but 53% occurred during the 
June - August period, with an additional 16% recorded 


during the "shoulder" months of May and September. Monthly 
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FIG. 4.1: MONTHLY TEMPERATURE AND 
TaGGe) PRECIPITATION 
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values varied considerably from year to year, especially 
during summer. February, the driest month, averaged 16.4 
mm, and recorded extreme values of 2.6 mm and 34.3 mm 
Extreme values for the wettest month, June, were 23.5 mm 
and 176.4 mm Longer records are therefore required be- 
fore true monthly means can be obtained, but the data are 
sufficient to serve as rough approximations with which 


data for 1976 may be compared. 


4.2.3 Sunshine 


Mean monthly sunshine hours varied from 72.1 hours 
in December to 310.7 hours in July. The most striking 
feature of the sunshine hours distribution curve (Fig. 4.2) 
was the comparatively low value for June, which, as the 
month of greatest daylength, would be expected to exceed 
values for July, May and August. The cause for this de- 
cline becomes readily apparent upon examining the precipi- 
tation distribution. High precipitation values also appear 
to affect September's sunshine input, when compared to 
October's; October's values were higher 5 out of 12 years. 

As with the aforementioned parameters, mean sun- 
shine values varied considerably from year to year. Ex- 
treme values for December were 50.2 and 92.3 hours, while 
for July, they were 247.1 and 348.7 hours. However, June's 


position with respect to July was fairly consistent, since 
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it had the lower value 10 out of 12 years. For May and 
August, there was more variability, as June exceeded their 
totals 5 and 4% times respectively. Although there is some 
variability in mean monthly sunshine hours, the rankings 
of the months are more consistent (May and October are 
exceptions), and make a reasonably valid basis for com- 


parison with 1976 values. 


4.2.4 Soil Temperatures 


20 cm soil temperatures at G plot were warmest 
overall, especially during summer (Fig. 4.3), with an 
annualemean of 6.2°C (lables@i sles ic Coldest at 5.1°C, 
due to cooler conditions in.-spring “and fall. 5B and F 
were much colder than the other plots during winter, with 
B recording a temperature of -5.9°C in January. It is 
interesting to notice one common feature for the plots, 
namely that during the transition periods between hot 
and cold seasons (spring and fall), the plots reacted 
almost homogeneously, while during periods of consistent 
heating and cooling (summer and winter), their monthly 
mean temperatures diverged. For example, mean monthly 
temperature differences of 5.4°C and 3.0°C occurred in 
January and July respectively, but during May and October, 
the ranges of values between the plots were only 1.5°C 


and 1.2°C respectively. This indicates both the influence 
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FIG. 4.3: MEAN 20 cm SOIL TEMPERATURES 
JAN.1970-JUNE 1973 
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of changing general weather conditions on the micro-climate 
or each plot, and the ability of each surface cover to 
modify its own micro-climate during periods of more stable 
weather. In other words, these covers can exert some in- 
fluence on their 20 cm. soil temperatures if general 
weather conditions remain relatively stable for a pro- 
longed time period, such as 2 - 3 months. 

These mean values, however, were calculated from 
only 3 1/2 years of records, and cannot be considered long 
term means. It is more difficult to examine soil tempera- 
tures than air temperatures over a short period because 
of the influence of general weather conditions, overlying 
cover, and the soil itself, especially its moisture level, 
which may not be the same from one year to the next. How- 
ever, as with other parameters, the months maintained 
similar ranks during the entire period. The exceptions 
are July vs. August and December vs. January. Ihese are 
respectively the periods of maximum and minimum temperatures. 
This is an indication thay the time of occurrence of annual 
extreme values is not consistent. 

At 100 cm, the influence of general weather and 
surface conditions on soil temperatures is somewhat re- 
duced, and the seasonal changes are less pronounced than 
at the shallower depths. Mean values below freezing oc- 
curred only on F, B, and PS plots (Fig. 4.4), while @ and 


S plots recorded brief periods of below 0O°C temperatures. 
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FIG. 4.4: MEAN 100 cm SOIL TEMPERATURES 
JANMP7OIUNE 1973 
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W has yet to record mace ee temperatures at this depth. 
As with the 20 cm curve, the 100 cm curve showed some 
divergence of values during winter and summer. However, 
major variations could also be seen in May, while April 
appeared to show a reduction in temperature range. Annual 
meansevaried from Seo cGestor. 6s  toso, 0°70. tor c) (lables... ). 
resulting from much warmer monthly values for G during 
late summer and much of autumn. 

The time of occurrence of annual maxima and minima 
at the 100 cm depth was more consistent than at shallower 
depths. At Hllerslie, the months were August and February, 
respectively, with the exception of W and G plots, where 
the minimum occurred in March. In fact, monthly ranks 
were more consistent here than at 20 em for all 6 plots. 
This curve also reflects the time lag that commonly exists 
between surface weather events and their resulting effects 
at depth. 

One may note that mean annual soil temperatures 
exceeded mean annual air temperatures by more than 4°C, 

It is apparent that the cause of this divergence was the 
winter situation, and the insulating qualities of the 
overlying snow and vegetation cover on the soil. Means 
calculated for G (used since air temperatures were re- 
corded over G) during the November - March and May - 
September periods are shown in Table 4.2, along with 


mean air temperatures for these periods. The difference 
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in winter temperatures was greater than 10°C at these 
depths, andeit Ve tantiicaupated’ that similar conditions 
exist within the tree plots as well. 

Although it will be difficult to make valid com- 
parisons between the above means and 1975 - 1976 values, 
the monthly rankings and time of minimum and maximum 


temperatures can be used as guides for later commentary. 


iw2.5. Fan Evaporation 


A 6 year mean was calculated from data for 1969, 
1971-673, and, 19754 76.data. There are data missing 
fors May, 1972, and 1975, June, 1975, and September, 1971 
and 1972, and no readings were ever taken during April 
or October. At best, these data provide a very rough 
first approximation, rather than a usable mean, so only 
a cursory examination will be given. 

May recorded the highest monthly mean value, 
6.39 mm/day (Fig. 4.5). Steadily declining values 
occurred during the remainder of the summer. This could 
be indicative of more humid conditions during summer, while 
September would show lower values because of its cooler 


temperatures. 
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FIG. 4.5: EVAPORATION 
CLASS ‘A’ PAN 
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4.2.6 Relative Humidity 


Weekly mean values for relative humidity are listed 
in Table 4.3a. Regression values (Table 4.3b) for temperature 
are listed for comparative purposes and will be discussed 
later. Relative’ humidity svalues sforsPsSjnand PWare 
probably not very accurate, but they can be compared to 
one another on the basis of thermohygrograph calibrations 
(Appendix 1) effected after termination of the field pro- 
gram. Since correlations and standard errors are similar, 
it can be assumed that the calibrations reflect the instru- 
ment error with reasonable accuracy. Although the magni- 
tude of the error is uncertain, the direction (positive 
or negative) is consistent. 

The slopes and intercepts indicate that PS is the 
Most humid, and FP, the least humid plot. There) is a large 
spread in slope values, indicating that relative humidity 
differences between the three can be substantial. Since 
the temperature slopes do not vary in the same manner, 
temperature alone cannot account for these differences. 
Therefore, stand structure and density, and their effects 
on evapotranspiration, should be considered as important 
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of cold dry winters and mild moist summers, with rapid 
transition periods. High summer precipitation is further 
reflected in the low values for sunshine hours in June, 
and for pan evaporation in middle and late summer. 

Soil temperature records indicate that consistent 
differences between the plots have already developed, with 
values diverging in winter and summer, while surface covers 
appear to have little effect during the transition months 
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4,3 June, 19 a April, 2976 


The 1l month period preceding the start of the field 
program was warmer and drier than the mean conditions 
described above. Ten months had below mean precipitation, 
especially the September - November period and the months 
of January and March, all of which recorded values rarely 
exceeding 65% of the 12 year mean. It was the second 
driest September on record (4.4 mm) at the station, and 
the winter precipitation as a whole (November - March) 
was 20% below normal. 

Temperatures were close to normal during the summer, 
but were consistently above the mean from November - April. 
January and April both recorded new monthly absolute maxi- 
mum values. Sunshine hours were close to normal, except 


for November, which exceeded mean values by 50%. Evaporation 
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data were incomplete, but showed mean conditions for July, 
and below mean for August and September. Soil temperature 
extremes were recorded in July for all 6 plots at 20 cm 
depth (Fig. 4.6). However, as a result of the mild winter, 
minimum values were reached in December, except on the S 
plot which recorded its lowest monthly temperature in 
January. S plot was by far the coldest, and was the last 
one to thaw during spring. This was also true at 109 cm 
(Fig. 4.7). The difference between S and the warmest 
summer plots, G and F, exceeded 5°C at times, while during 
winter, S was more than 2°C colder than W at 100 cm, and 
5°C colder at 20 cn. 

100 cm depth temperature extremes were reached 
during August and March for most plots. Although W and 
B had their maxima in July, the time difference is not 
considered significant, since their values were only mar- 
ginally above August values. The colder conditions for 
March as opposed to February, were probably a consequence 
of the dry winter, which reduced the snow cover, and so 
allowed the winter cold to penetrate to greater depths. 

Mean values for the period reveal that G and 8S plots 
were warmest and coldest respectively (Table 4.4). This 
is quite similar to mean conditions. Although PS remained 
much colder than the open plots, PS plot's 100 cm mean 
temperatures were somewhat higher than S because of milder 


winter temperatures. Snow cover was probably a factor, 
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FIG. 4.6: 20 cm SOIL TEMPERATURES 
T(C?) JUNIE A77SDEP 1. 77S 
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SOURCE: SOIL SCIENCE DEPT., U. OF ALBERTA 
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FIG. 4.7: 100 cm SOIL TEMPERATURES 
HED) JUIN S127. 5°SEP lento 26 


SOURCE: SOIL SCIENCE DEPT., U. OF ALBERTA 
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since the open hardwood canopy would allow more snow to 
reach the surface than the dense closed canopy of spruce, 
especially during a dry winter. F and G recorded highest 
‘summer temperatures at both depths, while W was the warmest 
WIntereplot.. hic eicecimmlarn: tomprevious years. 

As the 1976 growing season approached, conditions 
were ripe for early planting. Soils were dry at the sur- 
face, Snow cover had disappeared early, and buds had al- 


ready formed on the poplars by April. 
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April recorded above mean monthly temperatures and 
sunshine hours, and normal precipitation values (Figs. 4.1 
and 4.2). B plot was covered with last year's stubble. F 
plot was invaded by stinkweed (Datura stramonium), and 
there were weeds on W as well. Leaves began to appear 
on PS, P, and W, while G was still a yellow colour. Soil 
samples collected April 28 - 29 and examined gravimetri- 
cally for soil moisture, revealed extremely dry surface 
conditions under PS, S, and F plots, while P and W were 
quite moist at this level (Table 4.5). At depth, F had 


higher moisture contents than allthe other plots, thereby 


lay mean values discussed herein refer to those de- 
scribed in Section 4.2. 
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providing an illustration of the usefullness of a fallow 
rotation for conserving soil water. PS had the lowest 
mecorded values, but if data for 5 plot had been avail- 
able below 0.5m, it would probably have shown similar 
conditions, since this was the case at shallower depths. 
The field program started in May, so data were 
available on air temperatures and soil moisture content 
for each plot. Temperatures recorded over G plot were 
used for the other open plots (W, B, and F), while the 
forested plots had their own screens. Open air tempera- 
tures were slightly above normal while sunshine hours 
were below normal. Weekly mean temperatures were cal- 
Gulated for Pyro, S, and G pilots commencing (on the 13th 
(Fig .“4.8, Table 4.3a), and they reveal a wide variation. 
PS was 2°C colder than G on the 19th, and the other plots 
were cooler by at least 0.5°C. leaf cover on P and PS 
plots was complete by the third week. B plot was 
cultivated during the second week, and was then seeded. 
By May 27, sprouts as high as 8 cm had appeared. F plot 
was tilled and weeded while G was still somewhat yellow, 
with some dandelions (Taraxocum Dens-leonis) present. 
Soil temperatures were below 0°C for S and F at 
100 cm in April (Fig. 4.7). By May, all soils were thawed, 
Dut.o was stillsrecistering aschurly l.4.G at l00ecm mand 
6.4°C at 20 cm. W and B were warmest at 100 cm, while B 


and F had high temperatures at the shallower depth. The 
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only unusual aspect about this period is B's 100 cm 
temperature being higher than that for G, especially 
Since their winter temperatures were similar. 

S0il moisture levels for 0.1m (Fig. 4.9) fluc- 
tuated drastically in response to rainfall and time of 
moisture measurement, but showed a general decrease during 
the first 3 weeks of May. This trend was reversed towards 
the end of the month due to increased rainfall. A general 
downward trend was apparent for all plots at 0.2 m, which 
is probably the best depth at which to examine week to 
week changes in moisture content (Fig. 4.10). Data for 
Goan (ric. 4.21) seemmrosindicate allage time: that is 
somewhat longer than at 0.2 m, but well within the time 
frame that this study covered. May values at this level 
showed a downward trend for P2, P3, W, and PS, while the 
other plots remained relatively stationary. These results 
are probably a consequence of the dry weather which occurred 
SisApril and early May. | 0.0m. angyl.Omm) (lize 4 legend 
4,13 respectively) generally showed little fluctuation, 
and were obviously not altered by short-term weather events. 

There was little change in the plots' ranks during 
the period. The 3 P plots and W showed highest moisture 
contents at 0.1m, with values ranging from 35 - 45 cm 
water rie soil (cm mt). S was definitely the driest at 


23 cm are. The order waS similar at 0.2 m, but G was 


wetter than PS, B, and F, and almost as wet as W. Moisture 


110 


-gurt (0 ge sihso ct acoved iailtudin EB08 
Lo mals te {eile af ourogeer ak vssaSlaa 


aninnh sascioet: Lagendy a bawore tod treme 


ehcewor borvever Gawiient etd? “vyan ‘to aioe seat.¢ 
(inteiss beewortal of sub danon itt Xo, tne 


[geoney A 
avinw a $-0 24 erolg tie qot Ineteqqs eaw bnett a 
of 4e5w of bre ot doite +e ndash teed oft eidetorg a 
wot etal °.(OL.% walt) tastnoo ervtetom mt wean a 
: pt ten? ont? gal 2 dtagthal oF goose (Li.# wont 


eae eff cicdtiw flew fud =m S 0 vs. naity seers 
herevoo vbeta vide tent: 


ef ofidw ,2i pe W .€% ,S4 10 pawrt Beiewqwob a. i na 
edivear gua.  .ytenolserve yievitsler hogisaes atolq Tadeo 


doliw tedd.mew CID aus to eanenpesnne s Yldadanq. ote 
sit Yoras bas cen: 


fevel @rae 32 gates vem 


ipevsivse 
bon Sie .eptt) 0-4 fee a 8,0 
; yoottautoult? elists bewora viteyones* Ayievisooqaet el, 
er nere Yotteew oiat<P tone. yw bohed ke ton yiavotvde oraw a 
oniewh eines *etelq ett at eyando efftli esw e1edt 
cqugeion tootghd Bowens W brs atolg FC oft Lpotveq ade 
es se bee mame um £.0 ta edmernod 


te teviah sdf yletintted ase & ; Ty mo) rhea +" 
Pela prigark appar ing ao 


iat 4 dae Xow 2A mad 


pie o>, aay 


ariel) 


(>) “g1D5ud [Sh \ 
\ 


je 


O'Or 


iVase Jags . r ; > 
ate yet o Tia) Niadatiial 


r a a 
at © ’ 
—_—— —_ he 
(| ree Y iF Jes oreren Me 


a 


LI2 


0°0 


Ol 


07 


Ore 


( WD) “dID3¥d 


(02) 6"F “OIY 


3 


0°0 


0'l 


07 


oe 
(ttt?) edloadd 


0°02 


0'0€ 


00” 


($02) 6'F “Old ce 


114 


00 
0°02 

01 

0% 
0'0€ 

O'€ 

(aU) dlosdd 
0'Or 


(4409) 6p “OH ie 


TS 


( W) “gID3¥d 


7 
if 


oo 
bes 
7 


116 


a 


0'€ 


((W>) “dID3ud 


(4405) Ol’ “D]4 


AINE INN 
el 8 € SO €S Ql el 8 ¢ 
i | 
q- 
Agate ae ee sm 
ane =. 2 Ne 
He ef —— xe 
i yy A oe me 
Ese to 3) ee 
Dae a 
~ A 
Ss See 
XX, ~/ 
~~ : 
~~) 
ic) 


AVW 


00! 


0°02 


“Ss \ 0€ 


0°Or 


d7 


LO) 1d3$ “ONV 
he O linet Soonn OSC cpa CeaiG ge Glau IZ 
00! 
01 
oz 
ayer href 0°02 
O€ 
( Wd) ‘d1D3ud 
0'0€ 
0'OY 
00S 
(44¥09) :OL'P ‘OI4 w 


118 


Oe 


(atd> ediosda 


(4499) :OL'p ‘OI4 


uD 


WS "0" FNLSIOW TOs :| 


Dies ole 


001 


0°02 


‘0€ 


0'0r 


— 
a ; 


a 


- 
—_ ee a 
a ee ee 
” - 


120 


0°0 


0’ 


07% 


Of 


tee eciosad 


Lier al 


Cy ices (rc ct 


cae 


=—_— 
= _-_— 


— 


AINE 


OO AS os 


00! 


0°02 


DOE 


~ 


C ‘Or 


0°0S 


(*$4O)- 11°F “Old 


~ | 
° 


—— 
wa" = 


— w 


ak 


0°0 
00 
om 
ic 
—1 0°02 
O€ 
(o> )edioada 
00€ 
0Or 
0°05 


(4409) FLY “Og eng 


Ze 


One 
( WD) *41D3¥d 


Gio ler Ol 


00! 


223 


Y 8°O-FENLSIOW TOS :zL“p ‘Old 


——— Lh 
—_— eA! F 7 > 12 
——_—--— a : 


. Me 
4 ° 
7 o-oo 
-—* 
mS 


——o 


124 


0°0 


0'l 


07% 


Oe 


( WD) “gIDIad 


(4U0>) Zi" “Old By oe 


125 


o€ 
ey DERE 


(HOS) Sly “O14 


Ph. 


pion 
Wee) Fit Patten see a 


a. ae 


126 


Ot 
(°W>) "diD3ud 


(4492) ZL"y ‘Old 


N27, 


O€ 


( WD) “dID3ud g 


Se -_—=— 


— —_— 


aOR) USIOVVEl @Se Cle pale iso oy aa 


ae vo ee oe 2 

; a fa sa ——< _ 

vo : a - BALTIC ‘ 
> ; . a iv ey 102 3 . 


ae ae "ay 


128 


0'0 


0? 


O€ 
( W>) “dID3ud 


CEOs 2 Se 


i ee 
a 
> 
o Si Fee 
‘£7 = 
= 
< ie WES 
~ - - 
oo i yeas? 
* = - 
_ 
~ = a” %%; 
Pal 


rom a ae 


atid BAG Gy ae 


4. > 7 


- 
. 
i 
7 

i 

4 


contents were higher for most plots at this depth, 5S and 
F being the exceptions. At On5 im, PS) wasinow one of the 
driest, though not vas dry as 5, while G and W were nearly 
equal. S, B, and P3 were wetter at this level than at 
Shallower depths, while all others were drier. PS and 

F were the only plots wetter at 0.8 m than at 0.5 mM, 
while all plots with data available at 1.0 m showed 
higher moisture contents. B was the driest pLOt at One m 
aidaleo me with Vales.off28ecm.imictand 30 om ioe Wee 
spectively, while PS and W were fairly dry at 0.8 ty SEF 


were as wet as the other plots at 1.0 m. 


4.4.2 June 


Above mean precipitation and below mean temperatures, 
evaporation, and sunshine hours were recorded for this month. 
The first two weeks were fairly dry, with most of the pre- 
cipitation falling in convective showers. The last two 
weeks were wet and cold, with daily precipitation exceeding 
15 mm on three occaSions. F saw a return of the stinkweed, 
as well as Canada thistle (Carduus arvensis) and other 
plant species. This condition persisted until the plot 
was weeded and tilled on the 15th. B and W also suffered 
weed invasions, B receiving small numbers of ball mustard 
plants (genus Brassica), and W getting hawk's beard (genus 
Hieracium). Moist conditions in P led to increased herba- 


ceous cover, and poplar fluff was filling the air as the 
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trees released their seeds. By the end of the month, B 
Sprouts had reached heights of up to 0.4 m. 

Air temperatures for S plot became warmer than 
the other plots and maintained an advantage of slightly 
less than 1°C over PS, the coolest plot. G and P were 
fairly equal during the period, and maintained positions 
between S and PS. While S had warmer air, however, its 
soil temperatures were Heard ye5scecoldemethanuer F, and 
W plots at both depths. S was relatively cold as well, 
more than 1°C colder than the open plots. 

Surface soil moisture showed general increases as 
a result of the late June rains. The P plots remained the 
wettest, while W was the wettest of the small DLOtS and 
pene dries tes At Ono me: depth, S became more than 10 cm 
m7 Grver thant, yo. andaw by the end: opetheimonth: 
Before the rains, PS had dropped to within 5 cm. we ae 
S, but greater throughfall here allowed PS to gain 5 em 
mt, while S was barely affected. At 0.5 m, G values 
Showed a sharp increase on July 1, which may indicate that 
this was the only plot where infiltration of the previous 
week's rains could reach this depth. Slight increases one 
week ieee can also be seen at P2, F, W, and S and these 
may indicate a slower and substantially lower infiltration 
rate. Considering the very small change in its 0.2m 


values, the increase at S is somewhat puzzling, but since 
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the increment was only 0.05 cm ee forwa Wietam lots. 
this little bump in the curve should not be considered 
Significant. At greater depths, a slight downward trend 


was apparent at B and W towards the end of the month. 


4.4.3 July 


Mean temperatures and sunshine hours, near mean 
evaporation, and slightly below mean precipitation occurred 
during July. Between the 2nd and the 28th, no days had 
rainfall greater than 5 mm, and much of the precipitation 
fell during local convective storms. Heavier thunderstorms 
during the last week brought an additional 20 mm of rain 
to the area. Weeds gradually returned to F, and by the 
15th, the plot was heavily covered with vigorous plant 
growth. its surface was dry and crusty, with cracks evi- 
dent throughout. The plot was tilled and weeded during 
the third week, but the dead plant material was left on 
the plot. Weeds had also returned to B, but were not as 
serious. The sprouts were still doing well, and by the 
end of the month had attained a maximum height of 0.9 mn. 
Heads had formed by the 15th and these were large, promising a 
good yield. G and W also showed healthy growth and only a 
slight weed infestation. 

Air temperatures for S were still the highest, 


exceeding G and P by 0.5°C and PS by more than 1°C during 
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the last 3 weeks. S however. continued to register the 
coldest soil temperatures at both depths. F and W were 
now about 2°C warmer than B, and G at 100 cm had also 
exceeded B's temperatures. PS was still cooler than the 
open plots, but maintained a 2°C advantage over S at 

100 cm and @a=) CGC eatecOmcms 

The dry weather produced a downward trend for 
Surface soil moisture on all plots, but there were some 
changes in rank. B was now as moist as some of the 
P sites while PS declined rapidly. Toward the end of the 
month, W, F, and G were almost as dry as PS, while S was 
by far the driest at 25 cm Sons SS ul ae lower than 
the open sites. The rains during the final week brought 
sharp increases to all the plots exXCep teem sonceo. 
These did not react until greater Drecipitation sre. 
during the first week of August. 

At 0.2 m, a steady downward trend was evident 
for all sites except S, which was already too dry to lose 
much more water under these temperature conditions. By 
the 28th, B was Nearl Vaasa rVyeaseoy recording a valueror 
20 cm aed PS and P3 were 6 cm aoe wetter, while the 
other plots were clustered around the 29 - SRE oh ae range. 
At 0.5m, S remained steady while B dropped far below 
25 cm oe and so became the driest at this depth. The 


growing barley was apparently beginning to tap the 


moisture reserves at lower levels by capillary action 
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extension of an active root system. This trend was 
detectable in small magnitudes at the 100 cm depth as 
well. 


yy. August 


Above mean temperature and precipitation, near 
mean sunshine duration, and below mean evaporation, were 
the major features of August's weather. Heavy rains oc- 
curred during the first week, with 26.5 mm talline sonmethe 
ist. The other days recorded precipitation higher than 
10 mm, but it was that first storm which seemed to have 
the greatest significance on soil moisture levels. Barley 
on B had started to turn colour during the last week of 
July, and was now ripening. However, by the third week, 
its condition had deteriorated. Weeds (mainly thistle) 
started to dominate, and the barley stand was suffering 
from the competition. Weeds had returned to F and W as 
well. 

The early rains were accompanied by high tempera- 
tures, and the weekly temperature curve shows an increase 
over the end of July. G and S were now about equal, while 
PS started the month 1°C cooler. The range decreased 
during the middle of the month, but by the end, their 
values had diverged again. G was warmest at this point, 


while PS was more than 1° C cooler, having recorded a mean 
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Lemperature or fa lo. Sait temperatures showed increases 
fOrea se! plots atepoth depths, with the exception of F at 
20 cm, which had a slight decline. Its Status as the 
second warmest plot behind G remained unchanged, as did 
the ranks of the other plots. S was still 2°C colder 
than PS and 4°C colder than G and F at the 100 cn. level. 
At 20 cm, S and PS were only 1°C apart, but @ and F were 
Still much warmer than both tree sites. 

Surface moisture levels increased Sharply during 
the first week because of the heavy rain. S now recorded 
values higher than F and W, and had the greatest weekly 
increase, which may be due to a lower surface evaporation 
rate than at the open sites. Values started to decline 
again, until the third and fourth weeks, when additional 
rain brought increases to P, F, PS. and B. Towards the 
end of the month, however, only the P plots reacted to 
the rain of the 26th. The smaller plots were declining 
and for 5 and G, the decline had been continuous since 
the first week. The 0.2 m level showed a similar trend, 
though S remained more than 5 - 20 cm m2 drier than the 
other plots. B showed no decline after the initial week's 
rains, which may indicate the lack of water use by the now 
ripening crop. P maintained moister conditions, but G, F, 
and W were now as wet as P2, and wetter than all the other 
sites. B continued to be the driest at 0.5 m, while the 


3 P plots held slight advantages over G and F. The wetter 
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plots showed a general downward trend while the dry plots 


remained steady, as did all the plots at lower depths. 


4.4.5 September - October 


The final five weeks of the observational program 
were warmer, sunnier, and wetter than mean conditions. 
Evaporation was slightly below the mean. The precipi 
tation figures are somewhat misleading, because of the 
uneven distribution of rainfall during this period. Only 
5% of the 44 mm recorded during September fell after the 
I3th, sand&rainfall during the first week of October was 
confined to the 3rd and 4th. Two days registered precipi- 
tation above 10 mm, and only 11 days out of 38 had measur- 
able rainfall. However, another source of moisture entered 
the picture at this time--dew, which was observed 13 mornings 
during September (Table 4.6), ana probably had some effect 
on surface moisture readings at the open plots. Dew does 
not regularly form in densely forested areas (Geiger, 1965), 
and indeed, that rule held true for Hllerslie. The sharp 
decline in temperatures, which is an important feature of 
September weather, led to higher values for relative humid- 
ity, as well as the frequent formation of dew, especially 
over long grass. Fungal growth increased as a result of 
the wet August weather, with moss growing in PS, and toad- 


stools in S. F was once more infested with weeds after age 
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TABLE 4.6 
DEW FORMATION IN SEPTEMBER 


24 hr. period 


hrs. ST Min (°C) Dew Pt. (°C) Dew (x) 
1 10 10 - 
2 10 9 : 
3 3 = = 
4 7 9 x 
5) 9 10 x 
6 7 6 - 
? 3 3 = 
8 2 4 a 
y 5 7 - 
10 6 ? a 
i bak 8 6 = 
12 6 5 ¥ 
13 4 8 x 
14 5 7 se 
15 5 9 x 
16 5 8 ‘i 
ee D 9 - 
18 5 9 x 
19 8 Zi - 
20 5 ( x 
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TABLE 4.6 (con't. ) 


24 hr. period 
ending 0900 


hrs. MST Ming (eC) Dewe etm Ch Dew (x) 
Fait 5 9 x 
22 6 8 x 
23 6 8 ne 
24 8 iu = 
25 8 9 as 
26 4 9 x 
tat 5 5 = 
28 2 8 x 
Zo 6 8 ia 
30 4 uh - 
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had been tilled in August, and dew may have formed on the 
plant cover, which could have facilitated a rise in sur- 
face soil moisture late in the month, when the plant 
cover would be more extensive. This will be discussed 
later. 

Leaves at W, PS, and P plots were turning colour, 
especially during the last two weeks) DUL leatetall snac 
not been completed before the field program was terminated 
Ong0ctober So. By this time, however, those leaves still 
on the branches were yellow or brown. The thistles on 
B plot were also turning brown as a result of some cold 
night conditions, when the temperature approached 0°C., 

G remained green, as did the weed cover on F. 


PS was still the coolest plot, while the other 


?) 


sites had fairly similar air temperatures until the 15th 
wienGo cooled to the same Jevelvas PS, 11,0 Gs) wae now 
the warmest, but by the end of the program, all tour plots 
had air temperatures within 0.5°C of each other. § had 
the coldest soil temperatures at both depths, 3°C below 
Gand F. Rankings remained the same for this period as 
for August. 

Surface moisture was affected by the first week's 
rains, as indicated by the increased levels during the 
first two weeks. An exception was F, which showed a 
continuous decrease until the third week, when the trend 


was reversed. G and P fluctuated quite strongly during 
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the month, while W and B were more Stabie, averaging 


and 40 em m2 respectively. S, PS, and F 


37 cm m_ 
were the driest while P, W, and B maintained levels up 

to 13 cm ee higher. It is difficult to quantitatively 
determine the effects of dew on these cEbey ela Sn ORS 
possible that the increase registered during the third 
week at F, could have resulted from frequent dew formation 
on the weeds during this week. Another plot which might 
have been affected by dew was G plot, which suffered a 
moisture decrease that was less than anticipated, con- 
sidering the week's light precipitation. Dew was observed 
at B as well, but the denser stubble cover could have pre- 
vented much of the dew from infiltrating into the soil 
before it could evaporate. 

Although there may have been some surface effects, 
the 0.2 m level showed a downward trend after NemLirst 
week's rains. B reacted very strongly to the ralnwe jn= 
dicating the reduction in water use by the dying plant 
cover. B was now the fourth wettest behind MiG, and, 
whereas 6 weeks earlier, only S was drier. S and PS con- 
tinued to record the lowest moisture revels SPvAt £0, 5am 
depth, B again showed a strong upward trend as a result 
of the late August and early September rains. By the 
second week, its moisture level had climbed from 22 cm om 
to nearly 30 cm ee surpassing both S and PS. The other 


plots maintained stable levels about 3 - 5 cm nay higher. 
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W had been somewhat drier during August, but showed an 
increase similar to B's, and was now at the same level 
ASpGy ely Mande ee eAt. OMOimemand® 180 m depths, there was 
little change during the period, although B showed a 
Slight decrease at 1.0 mn. This may be a consequence of 
earlier water use by the previously healthy crop during 


early and middle summer. 


4.4.6 Summary 
General weather conditions for Ellerslie in 1976 
(May - October) were slightly sunnier, warmer and wetter 
than mean conditions. Evaporation was Slightly lower. 
Seasonal change in vegetation characteristics was af- 
fected by a dry spring, which produced an early canopy 
at PS and P, and by a warm wet August, which may have in- 
fluenced the demise of B before harvest time.+ The warm 
dry conditions during September allowed PS to maintain 
its canopy until well after observations ceased. 
Significant differences in air and soil tempera- 
tures were readily apparent throughout the growing season. 
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ate Soil moisture combined with high air temperature 
may be another factor. As the soil dries, the actual 
transpiration rate of the plant should fall below the 
potential rate. Since the potential rate is dependent 
on temperature, the plants’ wilting point must also be 
related to temperature (Denmead and Shaw, 1962). There- 
fore, it is possible that the barley ceased to function 
before the August rains could influence plant growth, 
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Regression analysis (Table 4,3b) reveals that S DL0bas 
air temperatures were warmest during middle summer (when 
G plot registered a temperature of at least 15°C). Ps 
and P were consistently cooler during this period. At 
lower temperatures, the temperature difference between 
the plots declined. § continued to have coolest soil 
temperatures, while other plots maintained their ranks 
(with some variations) during most of the Summer. Soil 
moisture ranks were fairly consistent as well, except 
for B, which fluctuated according to its state of health. 
There were variations in response to rainfall, and this 
was due in large part to throughfall (e.g. S) and plant 
water use (e.g. B). Short term weather and plant growth 
effects did not influence the 0.8 m and 1.0 m levels, 
except during B's healthy growth period, which lasted 
until the second week of August. Frequent dew formation 
may have affected surface moisture levelspat F ana G, 
but it is difficult to determine its Significance. 

It is readily apparent that long term conditions 
greatly influenced soil temperature and soil moisture 
values, and the ranking of each plot in these two cate- 
gories. Air temperature was more dependent on short 


term weather. 
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CHAPTER V 
EXPERIMENTAL RESULTS AND ANALYSIS 


oat Introduction 


General climate and field conditions were examined 
in the previous Chapter setlere writ ic appropriate to con- 
Sider the results of the radiation and energy balance 
measurements and derivations for the various De Ovs. ean | 
bedo, global solar radiation, and net radiation are ex- 
amined first, with special attention being given the 
prediction of Rn from (Q+q) using eq. (2.6). Having estab- 
lished the quantity of Rn available for the sensible, 
latent, and soil heat fluxes, one may now discuss the 
variations of these fluxes. This is followed by an 
analysis of the Bowen Ratios. unergy balances are cal- 
culated based on weekly values derived from regression 
equations of Rn and Gs, and on calculations of LE. 
Measurement and statistical errors are then accounted for 
in an effort to determine the reliability of the energy 
balance calculations, and consequently, the validity of 
the possible micro- and macro-scale climate effects 


caused by future land use changes in Central Alberta. 


Dee Radiation Balance 
Sean Albedo and Global Solar Radiation 
aoe Yn viobal solar hadlation 


Daily mean values for albedo and (Q+q) are listed 
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in Table 5.1. Valves for June (Qtq) were slightly lower 
than for July despite a greater daylength. This is to be 
expected, since July had more hours of bright sunshine 
and lower rainfall. There was a sharp decline from July 
to August due to a decrease in bright sunshine hours, 
while September's values benefitted from a dry 3 week 
period with an abundance of clear days. 

Mean albedo values ranged from 0.065 for F plot 
in June, to 0.250 for G plot in September. S plot was 
the most consistent, while W varied from 0.196 in July 
to 0.104 in September. Since there was little change 
in S plot's canopy structure during the summer season, 
lack of albedo variation is no Surprise. Admittedly, 
the June data for W plot are based on just one 24-hour 
period # but July's values are only 1% higher than June's, 
and more than 9% above September's. The low September 
values for W, and PS, the other deciduous cover, were 
probably a consequence of leaf loss. 

B plot's values seemed to fluctuate with the state 
of the crop's growth. In June and July, the plants were 
young, and green in colour, so albedo values were high. 
Towards the end of Summer, however, the crop deterior- 
ated, its colour changing to a dull brown, and dark green 
weeds had infested the plot. Meanwhile, the healthy large 
barley field had reachea maturity, and its colour was 


golden brown. There were few weeds, and moist bare soil 
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was visible between the stalks, thereby reducing reflec- 
tivity to a value 2.5% below B plot in summer. 

F plot was affected by human activity. Higher 
albedo values at F indicated periods with substantial 
weed cover, and a dry crusty surface, while tillage and 
weeding resulted in albedo reduction. June albedo for 
G plot was affected by wet weather, which produced a low 
Val uerosfORLAtonthem2zth 

Overall? °GUplot andstie young barley crop had the 
highest reflectivity, 3 - 5% greater than W with leaves, 
and 8 - 10% above S, leafy PS, and dry weed-infested F, 
They were 10 - 15% above leafless PS and W, and moist, 
weed-free F. Considering the magnitude of daily global 
solar radiation, this variation in reflectivity could 


2 


amount to a difference of up to 40,000 Wm day + of Rn 


available for use by the various plots. 


Seoee Net Radiation--Diurnal Variations 
ao’ eto fonasviurnal Vyariations 


Diurnal variations of Rn for selected clear 24 
hour periods are illustrated in this section. These 
periods were chosen to maximize the. effect of the indi- 
vidual plots’ albedos on their respective radiation bal- 
ances, and to facilitate the detection of any measurement 
errors that may have occurred during the field study. Data 
for Gs and (Q+q) are included where available. Rn data for 


G plot are illustrated in all figures for comparative pur- 


ok cobet ice ee estate: cbedite: 
wol # vena poke pedeasy dows ‘iindaal ae 
gir aeons we | 99S odo $420 To auEEy: 
arid bad. qos yetsad peony: ad: bie Yolq ‘SoU LLasev0 : red 

mevact ‘tin W Godt Terao RE ~ © tivisosttbs seentgil © 
| bet pot ndehome? Ayah RAB RL Yaod 52 avoda SOL ~ 8 bas 
-betom ban: W-bnn- GF aeetiacd-ovetis Red + Df eapm at 
tadotp (sled td ~ebed Letgen" art grixebtenod- tT éertt boew ae 
biueo  YElviteetietand! mottaizey eias toltaIbaa ‘nf8 4 
pit te abe ate 9G O88" qu hovennd torte ‘* of ettoms 
hn ames - vetoiq: avolaay ont - siemiless eltdalisvs 


as aeelo weFohion a ak a teeele er 
quadt., emolsoon..citt, al bedextauls, orn, sholreq wor 
Fepecompesneneeses yee 


~fad enitaiig gvitosgeen | aledy mo wobed is. tatoiq ‘, we 
“or zat ea noirsetsh: ‘ead sist Ltiost of bMS -« isa 


gate perw300- evad yar | pale 3 1S 
ete seb eg pero sales ae 


146 


poses, since G was the base station. In addition, two of 
the figures contain data for Stony Plain (Q+q) and Rn. 
These will be discussed later. 

The calendar date given represents the 24 hour 
period ending at 0900 hours MST on that date. Therefore, 
the data from 0900 - 2359 MST were actually measured one 
calendar day before the given date. This was done so as 
to synchronize the small plot measurements with the daily 
meteorological 0900 MST observations at the Ellerslie 
Climate station. 

A note of caution is given before discussing these 
data. The graphs represent individual 24 hour periods; 
this is important when considering their interpretation. 
It should not be assumed that similar diurnal variations 
have been occurring for any specified time period before 
or after the given date. The purpose of this research 
project is to examine gross energy balance differences 
over several months. Detailed analysis for much shorter 
time periods requires greater instrumentation, and is 
beyond the scope of this study. 

Diurnal trends for PS plot are illustrated in 
Pigceeo UmAandes, oe mAcemichnt be expected, June 19 (Fig. 
5:1) featured much higher midday (Q+q) values than 
October 2 (Fig. 5.2). Upon comparing the two graphs, 
certain major differences in longwave emissions can be 
seen. On June 19, nighttime In values for G and PS were 


fairly similar, while in October, emissions from PS were 
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much greater, probably due to higher nighttime tempera- 
tures in the sparsely leafed tree canopy than in the 
open air. This could have caused a lowering of daytime 
Rn values for PS, so that its peak was now the same as G 
plot's, and its morning and late afternoon values were 
lower. Considering the fact that G plot's albedo was 
twice that of PS in October, and that on June #0, PS 


2 at 1300 MST, this 


exceeded G plot by more than 200 Wm 
change in rank is remarkable. Gs Values sat Giplotrex= 
ceeded values at PS during both days, while there were 
lower nighttime emissions at PS. Both occurrences il- 
lustrate the influence of the deciduous canopy in shading 
the ground during the day, and returning longwave energy 
to the soil surface at night. 


S plot's diurnal trends for June 30 and September 


26 are illustrated in Figs. 5.3 and 5.4 respectively. The 


results for June are similar to those of Fig. 5.1, although 


the midday Rn difference is smaller, and S$ ikea ch ie Yer ein 
time emissions are consistently greater than G Dios. 
However, the September situation is radically different 
from Fig. 5.2. Rn at S plot remains higher, and its Ln 
emissions are lower than G plot's. Considering the leaf 


loss on PS plot and the dense canopy of S plot, one may 


conclude that differences in canopy structure are probably 
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the major factor in this situation. 

Rn slightly exceeds (Q+q) during the late after- 
noon of September 26. This could be due to a short period 
of positive In, or measurement error. Error possibilities 
are discussed in Section 5.2.5. 

W plot data for May 31 and July 24 (Figs. 5.5 and 
5-6 respectively) reveal that In emissions are similar to 
those of G plot. However, Rn values are higher during all 
daylight hours (as opposed to PS and S$), thereby illus- 
trating canopy structure differences. Its rough exposed 
surface allows greater absorption of low sun angle radi- 
ation than the relatively smooth G plot. Diurnal varia- 
tions of Gs are similar to variations on PS Dio t. 

B plot data for May 27, July 4, and September 17 
(Figs. 5.7 - 5.9 respectively) show seasonal variations 
in Ln emissions, when compared to G plot. On May 27, Ln 
emissions were greater at B. On July 4, the difference 
between B and G was smaller, and on September 17, the two 
plots had similar In. Gs changed as well, with higher 
values being recorded for G plot in May and July, and for 
B plot in September. This illustrates the effect of 
changing plant cover and water use on B plot. In late May, 
surface moisture started to decline until heavy rain fell 
in late June. Heavy rain in late summer created wet sur- 
face conditions since plant water use had decreased. Moist 
soil has a higher thermal conductivity (K) than dry soil 


(Sellers, 1968, p. 128). According to the Gs equation: 
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Gs = -K a (5.1) where = = temperature 


gradient 


this should result in an increase in Gs in September since 


= In July and late May, 


heat storage capacity is increased. 
Gs values at B plot were reduced by the overlying plant 
cover, and low moisture levels, respectively. Nighttime 
soil heat losses were reduced somewhat due to the overlying 
plant cover. 

F plot data for June 2, July 18, and September 18 
(Figs. 5.10 - 5.12 respectively) illustrate the effect of 
F plot's variable weed cover on albedo, and consequently 
its Rn values. On June 2, maximum Rn values were as much 


2 


as 150 Wm“ higher than G plot's, while on September 18, 


the difference was only 50 wn”. in emissions were gener- 
ally greater at F plot, although on July 18, the heavy weed 
cover had an influence in equalizing the nighttime Ln on 
the two plots. 

Certain other characteristics of the above diurnal 
variations can be noted. Many of the Rn curves show a 
peak at night. Time of occurrence varies with the plot, 
and with season. The seasonal variation is probably a 
consequence of prevailing weather, but the plot variations 
are interesting. F, PS and W recorded peaks 1 - 2 hours 
eavermetian G plots (higs.e5. 10,95, leand 5.0) cumuuneme 
June 19, and July 24 respectively. B and G recorded 


lthis is similar to energy transfer in a calm water body. 
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maximum Rn at similar times (Figs. 5.8 and Sis Csvren lens 
showed little temporal variation with G. There does not 
appear to be a connection with nighttime Gs values, so 
this could be illustrating micro-scale differences in 
relative humidity. These differences are discussed in 
Section 5.2.4. Possible errors are examined in Section 
5:2.5, including potential instrument calibration errors. 
Stony Plain (Q+q) and Rn data were inserted in Pigs 


and 5.8 for that reason. 


SYN aris. Net Radiation--Weekly Totals 


Net radiation for G plot, the base station, was 
calculated from direct measurement during May 21 - August 
3 and September 15 - October 7. During the intervening 
period from August 4 - September 14, Rn was derived from 
(Q+q) using eq. (2.6). Slope and intercept values are 
listed in Table 5.2, along with correlation (r), variance 
Gay" standard error (S.E.), and level of significance 
(SIG). Slope and intercept fluctuations are in large 
part a result of variations in longwave emissions at night, 
and daytime albedo. Correlations and variance are very 
high, however, as are significance levels. Standard 
errors are less than 10% of average hourly values of 
(Qtq). A reciprocal regression was computed as a check 


against unusual variations in the relationship between 
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Rnwands(Qtq).— If ithe relationship is genuine, the slope 
should approximately equal the inverse of the slope 
using eq. (2.6). The inverse is listed in the table, and 
the slopes are within 0.06 of agreement. 

A variation of eq. (2.6) was employed to predict 


Rn for the rotated plots, based on Rn data from G plot: 
Rn (rotated) = a{Rn(G plot)} +b 552) 


Regression slopes and intercepts were computed at various 
times during the study period from paired data which con- 
sisted of simultaneous hourly measurements of Rn at G 
plot and at the rotated station. These were interpolated 
during intervening periods when the rotated set of instru- 
ments was measuring another plot. 

Regression slopes and intercepts to be used in 
eq. (5.2) are listed in Table 5.3. With one exception, 
B plot on June 22 - 24, all correlations are greater than 
0.90, and at least 80% of the variance is explained. 
The majority of standard errors range from 20.0 - 50.0 


B and one above 


Wm“, with only four higher than 60.0 Wm7 
80.0 Wm7*, Ona percentage basis, most of the standard 
errors are similar to those calculated for eq. (2.6). 
An instrument malfunction during June 22 - 24 was the 
cause of the poor correlation and high standard errors 


calculated for this period. The slope and intercept were 


not used in subsequent calculations. 
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Slopes ranged from 1.10 - 1.48. B plot registered 
the lowest ones, an indication that its albedo was the 
closest to G plot of all rotated stations. PS plot had 
the highest slopes, as well as high negative intercepts 
which slightly reduced Rn values. However, Rn values 
remained higher than those on S, W, and F plots since 
intercepts accounted for only 0.3 - 0.5% of Rn. Results 
for F plot were influenced by the same factors which 
affected its albedo; a dry crusty surface and extensive 
weed cover, both of which established themselves rapidly 
while the plot was undisturbed. Continuous tillage and 
weeding would probably have reduced the albedo, resulting 
in higher Rn values, and hence, higher regression slopes. 
Lower albedo values for W plot during September affected 
its slope, while S plot's slope was influenced more by 
changes in G plot's albedo rather than its own. PS plot 
represents a more complicated situation, due to the nature 
of its stand structure, and seasonal variations in leaf 
canopy density. Multiple reflection and absorption al- 
ways reduce the albedo of a deep stand of vegetation. 
However, the albedo of an individual leaf is much higher, 
so leaf angle, sun angle, and foliage arrangement can 
temper the albedo reduction. 


Calculated mean daily Rn for 5 - 9 day periods, 
using eq. (5.2), are listed in Table 5.4 and illustrated 


in Fig. 5.13. F plot had the highest values in late 
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May and early June, but PS plot took the top rank on June 


10, and maintained this position throughout much of the 


season. G plot had the lowest values 


period. The difference between G and 


low of 12,000 Wms” : 


during the entire 


PS ranged from a 


day ~ at the beginning of October to 


60,000 Wm” “aay7? during the first week of September. 


Significant differences occurred during much of August, 


as PS absorbed 20,000 - 40,000 Wm~2 


day 


= more energy than 


W, F, and S plots, which were themselves 20,000. -— 30,000 


wm” “day7? higher than B and G plots. 

which must have had an influence were 
surface conditions which affected the 
W plot ‘s@rankimeaduringsthe.firstetwo 
period illustrates the potential this 


sorbing radiant energy if weeding and 


Again, the factors 
plant structure and 
plots' albedo values. 
weeks of the study 
surface has for ab- 


tillage are main- 


tained. Although this increase in energy might have in- 


creased evaporation, the weeds could very easily consume 


as much if not more soil water, and their roots would 


be able to draw water from greater depths than would evap- 


oration alone over a short time period. At this point, 


the idea is purely speculative. However, a detailed 


examination of Rn differences between 


the plots, and 


the partitioning of Rn into its various components, 


especially LE, will go a long way towards settling 


this and other issues. 
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See Radiation Balance 


Radiation balance was calculated for all plots*on 
the basis of albedo, global solar Vadtacion seandanet 
radiation data presented in previous sections. Results 
acemuisted inet abliemsn 5. 

Differences in the prevailing weather of each 
selected period are apparent in the resuits. In values 
were close to 0 in August, illustrating the influence 
of this period's humid conditions. Comparatively drier 
conditions during the other 3 periods allowed greater 
emissions of longwave radiation from the various sur- 
faces. 

Significant variation between the plots is ap- 
parent as well. G and F plots emitted the greatest amount 
of longwave radiation in June and August. They had to 
share top rank in July with S plot, while most plots had 
similar longwave emissions in September. The plots with 
the lowest emissions were B and PS in June, W and B in 
July, and PS in August and September. 

Upon closer examination of the calculations, it 
is readily apparent that the main factors which are in- 
fluencing these emission differences are stand structure 
and growth characteristics. PS plot's canopy, which al- 
lowed greater radiation absorption to take place, is also 


reducing longwave emissions. The various leaf layers 
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readily aneann these emissions and reradiate them to the 
cooler surface below. This effectively reduces vertical 
temperature gradients above the canopy, while a tempera- 
ture inversion develops below, especially during clear 
daytime periods. The positive value calculated for August 
is probably a result of the combination of prevailing 
humid weather and the above structural effects.) @HiadWair 
temperature data been available for the canopy layer it- 
self (rather than the plot's floor), temperatures warmer 
than the other plots would have been recorded this month. 
Previous studies in larger deciduous forests have shown 
the crown to be as much as 6°C warmer than the floor, de- 
pending on time of day, species, and stand age and density 
(Geiger, 1965). Although this large magnitude of tempera- 
ture gradient would not happen in PS plot (since the stand 
is small and young), the radiation balance results indi- 
cate that if the plot was a larger stand, substantial 
temperature increases in the crown area would result. 
This is definitely a topic for further study at Ellerslie. 
Similarly, S plot's canopy structure is influencing 
its longwave emissions and subsequent air temperature. How- 
ever, Stand height is lower, and density is greater than at 
PS plot. The air temperatures here were actually measured 
in the lower part of the crown. Consequently, these were 
consistently higher than PS plot's temperatures. The 


greater longwave emissions are an indication of the 
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comparative dryness of the air above S$ plot, as opposed 
Om om enon omnia t surprising considering the fact that 
PS is probably transpiring at a greater rate than Ss. 
Soil moisture utilization is illustrated in Section 

BAS hae 

B plot's changing state of health is well illus- 
trated in these results. In June and July, the actively 
transpiring crop had one of the lowest emission rates. 

By September, longwave emission was similar to the other 
plots, indicating a sharp reduction in water use. 

W plot, the other deciduous cover, had lower 
emissions than PS plot during July, but ranks were reversed 
in June, August and September. This indicates a difference 
in length of active growth period. While PS remained active 
until well into the autumn, W started to lose its leaves 
during the first week of September, and its soil moisture 
levels showed a large increase during the last week of 
August. 

G and F plots showed similar emissions throughout 
the study period. G plot's emissions during July and 
August were higher due to F plot's heavy weed cover and 
lower albedo, which combined to produce higher water 
losses due to transpiration and evaporation. The weed 
cover also reduced F plot's emissivity. At wavelengths 
shorter than 2.4 microns, many plants have lower emissiv- 


ities than moist bare ground (Sellers, 1965, p. 41). 
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SILT 5 Errors and Analysis 


Several sources of potential error exist when 
measuring and deriving radiation balance components. 
These include instrument placement, recorder reference 
points, calibrations, data extraction, and calculation 
errors. 

Accurate data require that the instrumentation 
be properly maintained, and that the radiation sensor 
should have an unobstructed field of view. The height 
of the sensor varies with the size of the DLOt, sande. 
the plot is small, this height must be reduced so that 
the instrument senses only the desired area. The only 
plot which posed certain difficulties in instrument 
placement was PS, since the instrument tower did not have 
sufficient height to completely clear the canopy. The 
instruments generally had unobstructed views of the sun 
during all periods but those with very low sun angles. 
Since radiation intensity during low sun angle periods 
is small compared to total incoming radiation, this 
error is not considered significant when calculating 
daily and weekly totals. In fact, a reduction in in- 
coming radiation to the sensor should have decreased 
Rn values, and results show that PS had the highest Rn 
of the six plots, for most of the study period. The 


October diurnal graph (Fig. 5.2) shows Rn at PS to be 
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lower than at G. However, Table 5.3 reveals that the 
Slope of the regression line is still greater than 1.0. 
As previously mentioned, the diurnal variations should 
not be considered as representative of much longer time 
periods, especially in autumn when PS is undergoing rapid 
leaf loss. 

Albedo measurements at F plot posed a different 
problem, since F was affected by human activity such as 
weeding and tillage. As explained in Section 5.2.1, this 
factor is important to the plot's albedo, thereby making 
the date of measurement equally important. This will 
anGeCimtieeresultssofseq. (2.4), specifically In, putathe 
exact magnitude is unknown. A measurement error greater 
than +5% is not expected. Thus, the maximum potential 
error would be found in June's Ln values, since (Qtq) - Rn 
is greatest. However, the seasonal total of this error 
should be small, since measurements were effected on 
widely separated dates (see Table 5.3) and over various 
surface conditions. 

All instruments were factory or laboratory cali- 
brated, and these should have been fairly consistent 
throughout the measurement period. In fact, the Eppley 
pyranometer measuring albedo was brand new. The digital 
counter being utilized to monitor (Q+tq) was abandoned in 
favour of data extraction from the strip charts when it 


was Suspected that the digital was not functioning 
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properly on days with low incoming (Qtq)e me iheedtaital 
error is illustrated in Appendix 2. 

Potential calibration errors were checked by 
comparing (Q+tq) and Rn over G plot with published data 
obtained from Stony Plain (lat. 53° 33'N, donee Tt 06° wW)* 
While planning the field program, the idea of using Stony 
Plain data directly was discussed, but discarded in favour 
of direct measurement at Ellerslie because of potential 
differences in prevailing weather and its associated cloud 
cover. Table 5.6 illustrates the data difference between 
the two stations. The Rn differences could be indicating a 
combination of variations in prevailing weather and station 
surface albedo and Ln, as well as instrument calibration 
error. The surface effects would not be a factor in the 
(Qtq) differences. The significance of the error on this 
difference is hard to determine. It is quite likely that 
theregis a ‘calibzation errom, bat it would gappearmthas 
this is within acceptable limits, considering the fact 
that instrument errors of 5 -.15% are common. Subsequent 
detailed analysis will employ 2 week totals so that random 
variation can be reduced, and since the same data are em- 
ployed for all rotated plots, this will not affect any 
differences calculated between them. Interestingly enough, 
if the Ellerslie (Q+tq) were to be raised by 11%, albedos 
would be lowered 2 - 3% since the denominator of the 


fraction would be increased, and this change could create 
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somewhat questionable results, especially for PS and G 
plots: 

The Rustrak strip chart recorders also presented 
a problem. Their calibrations and reference points were 
periodically checked while the recorders were operating 
in the field, but since these are battery powered, the 
reference points may vary with the strength of the bat- 
tery. Battery power was continuously monitored, however, 
and recorder reference points were checked so as to ensure 
proper strip chart analysis. 

Other errors in strip chart analysis are human 
in origin. Random error is expected, but if it is ran- 
dom, the sum total of the error over daily or weekly 
periods should approximately equal 0. Comparisons of 
strip chart analysis for the same 9 day period, completed 
on separate occasions by the author and Dr. Hage of the 
Geography Department, reveal a difference of 0.118 chart 
units on a chart scale of 100 units. The error is thus 
considered to be random. 

Calculation errors are also human in origin, but 
these should be random as well. Major errors would re- 
veal themselves upon comparison of measurement results 
against known norms. 

Having considered the error possibilities, one may 
explain the radiation balance data for each plot, at least 


on a comparative basis, with those data obtained from other 
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plots. Equal amounts of global solar radiation were 
available for each surface cover at the same time. But 
Table 5.5 reveals marked differences as to how much of 

this energy was stored for later use by the various sur- 
faces. PS plot consistently recorded the highest Rn values 
and the lowest outgoing longwave radiation emissions except 
for the first and last two weeks of the SLUudyr. 2 and G 


plots had the highest emissions. Rn differences were 


2 7 


quite large, and in some cases exceeded 60,000 Wm “day. 
An imbalance has now been created, and this could affect 
water use and the other components of the energy balance. 
Therefore, it 1s anticipated that major “differences in 
energy balance will be calculated for the various plots, 
and these may contribute to the differences in air and 


soil temperature and soil moisture which were observed 


during the study. 


bao Energy Balance 


Sy. Shel LE 


APW and cumulative 4PW are illustrated in Figs. 
5.14 and 5.15 respectively. Results of eq. (2.12) are 
Pistedein Table 5.7 and illustrated in Figuss. 16.) In- 
cluded in the figures are weekly rainfall data, and in 


Fig. 5.16, dates of human activity on F plot. 
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vl general, all plots had reasonably similar LE 
values at the commencement of data collection. Most plots 
reached a peak at the end of June. However, subsequent 
calculations show the plot results diverging as the sum- 
mer progressed. Although large scale divergence is not 
readily apparent in the APW curves (except in mid-July), 
the cumulative curves reveal large variations between the 
plots at the start of the study. At this point, it is 
interesting to note the similarity in trends that most 
of the curves appear to follow after June 17, although 
it must be remembered that the plots started the growing 
season with different amounts of soil moisture, which 
would affect both APW and cumulative APW. However, the 
Same cannot be said for the LE curves. Reactions to 
weather events were inconsistent, and overall trends are 
difficult to see, except for a general downward trend 
from September 2 until October. Similar fluctuations 
exist in the LE/Rn curves (Fig? 5.17). 

Regression analysis revealed few significant 
linear relationships between the various plots’ LE and LE/Rn 
values (Table 5.8). Correlations were fair for PS vs S$ and 
PS vs W, but it would be difficult to arrive at any conclu- 
Sions based on these results, especially when one considers 
that the magnitude of the standard errors in the LE re- 


lationships often exceed 30% of the 20 week mean, and 
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less than 50% of the variance is explained for all but 
the PS vs S case. 

It is apparent that no Simple relationship exists 
between the plots' LE values, nor can LE be linearly re- 
lated to the plots' respective Rn values. Thishisran 
agreement with Davies (1972), who states that LE/Rn is 
not a general indicator of soil moisture. Therefore, 
more detailed analysis is required before any conclusions 


regarding soil water use can be reached. 


epee Gs 


Soil heat flux for G plot was calculated from 
direct measurement during May 21 - August 12 and September 
17 - October 7. During the intervening period from August 
13 - September 17, Gs at G plot was estimated by employing 
Gs/Rn ratios calculated before and after this period.? 
Data calculations followed a procedure similar to eq. 


(5.2) using regression slopes and intercepts: 
Gs (rotated) = a{Gs(G plot)} +b (5830 


Although daily totals had to be estimated during the inter- 
vening period, hourly data were available at certain times, 


and regression slopes were calculated from these. 


Lineage ratios are listed in a later section. 
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Regression slopes and intercepts are listed in 
Table 5.9. Correlations are generally lower than those 
listed in Table 5.3, but most of them are acceptable. 
Several individual periods, such as PS on September 3 - 4 
and F on August 29 - September 2 were hampered by recorder 
problems, so their sample sizes were small, consequently 
affecting significance levels. S plot had two correlation 
values below 0.8, but these were probably due to physical 
factors rather than instrument error. The coniferous 
canopy's influence should vary only Slightly, while Gs 
values at G plot fluctuate. This variance is magnified 
because the magnitude of Gs is so small. Therefore, large 
prediction errors are not anticipated, but to improve ac- 
curacy, the results for July 27 - 30 will not be used in 
subsequent calculations. 

F plot, as expected, had the highest slopes because 
of greater soil exposure. 5S plot had the lowest ones for 
the opposite reason. Most of the intercepts are within 


© of 0.0 Wm 


5 Wm , indicating that there are no major 
temporal variations between the plots diurnal trends. W 
and PS plots showed large slope increases towards the end 
of the study period, and this was probably due to leaf loss, 
which increased surface exposure. F plot's slopes were 
probably affected by weed cover density and soil wetness 
The latter point’ Is more important for F plot than for 


any other cover, since there is little or no overlying 
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vegetation that can intercept precipitation and reduce 
direct throughfall. B plot was also affected by surface 
wetness, but this varied with plant density and water use 
rather than precipitation frequency. During the summer, 
the healthy plant cover shaded the surface. Towards 
autumn, much of the plant cover was removed, thereby 
exposing the surface to solar radiation and wind. These 
events must have affected soil thermal conductivity at 
both B and F plots. Further discussion on this topic 
will continue in a later séction. 

Soil heat flux data are “lis*ted™in “Tabte 5£10¥and 
illustrated in Fig. 5.18. While its magnitude is quite 
small compared to the other components of the energy 
balance, significant Gs differences between the plots 
may, in large part, account for differences in their re- 
spective soil temperatures, which are illustrated in Figs. 
4.6 and 4.7. 

Upon examining Fig. 5.18, one notices certain 
consistent patterns. F plot's values are the highest 
(except for one week in June) and G plot's are the second 
highest from May 27 - September 16. At that time, sig- 


nificant differences disappear, and the daily soil “energy 


2 


balance" hovers fairly close to 0.0 Wm ~ as the winter 


season approaches. Another interesting point is that in 
May and June, differences can be discerned between the 


remaining open plots (W and B) and the two wooded plots. 


coxtava ys Hergerts wise lean telg @ <h | 

eay Tetaw bas ytteneb sea Asn botany abe da are 

ebiswoT aa - 

ydetent ,bevomes ssw TEVOD » dng efit to doum ameyue 

evetT  .baiw bra nolaethes, safoy of gastave (aut gnizogxs — 

s Wivisoubmos Lanrvestt live bet oatta oved teum atnsve 

siqot eidt so roleswoaib rode arolg & bos @ isod 

— wold oes setel s ai evnitnoo Iitw 
brie on. eids? ai befell aus stab xult teed Lio® 

etivp ei sbutingem afi ofiniW 8L.2 «att nt betasteutil 

ynvese of to atnerogmeo ‘redto edt. ot bexsgmon net 

atoig att ‘neow?ed seoneteltib 2D tasoltingte 

-o1 tnd? gl esomerettlb Jol fnlones «778¢ omnel ab am 


ealt et bazerteatth, run tat jeoansattsrmes ‘Lhoe ev 
: id ot ba Oat 


~« 


_ ee ‘ — ial _ a 


_ 


aoa _ 


nietzeo weolten gn0-45Lse eltushaimxe mony 
$ oediyhet) ort? ous ‘aouiav s'dolq remnettag snohelenes 
brrooes outs ou ately 5 bre (enul mi seew ono Zot taeoxe) 
-pin walt tang th 994 pedmerqee = $5 vet nowt teedatd 
yyreas” five yileb orf) bne weeqqnaib seonoaettin snssitta © 
i ee seo ls ‘frist asoved “eons led 
ni tart ai tnhog gititessetal serttonA I canta tia’ / 


oe eee eel 


owt eft bas (@ bie W) atolq: 


201 


TOC 2S 0°s q°€ 
9°22 0°9 GS 9°47 
E'St O°t 4° € U's 
4S¢ (hers Ona CaS 
Soy, 8h 7S eg 
Q* ce 9°6 6° UT 6°S 
5 er 6 2 OL oe 
EL hers €'g 4°€ 
Tec 9 u°9 Lene 
S°6T G°6 4° 6 ge 
OL 6°S T°9 ge 
Ge GSI GST 6°4 
d q M S 


(, Aep, UM OOOT) SD ATIVG NVAW 
OTS oSIEVE 


2°? 
ras 
6°¢ 
Cad, 
G<S 
2°6 
O°a 
O'S 
lees 
ee 
a2 


gS 


Sd 


€/9 
L2/S 
sutpug 


potseg 
a1eq 


ero) «SA ene 


’ . : a ae 7% ky ay 
pore 4 sh ah 4? 
% ee b 
Or ee . + 


wo: @ 


weve oa —nited on. 190s A. 


ye"r "0 


2s 


os 


boda’ kes weSOh oe hal a a aoe 


witnget 7 cunrh @ eet sy | | 
i Aagrt aan TOVee sents S80 Peas | je . 


it 
» releuern mrt re e ea al ea 2 
en one 2 OR sat | 
BS Go RBS aA BSE ey ra _ 

nae ah Ache unt Sins amc 


1% By? Leen? 
ao» tt 7 t+ 9 2 4a + : 
‘v etre fee ee ae nvieare oat 


ae ee a c 
estes nee trot tee 
> dn ow WD eo eee Oe | 

+ PEEERREt Salts 


se) 12? ey raete eae iii 


Rare a a a ae 
‘oa 


. Vip cele pein 


202 


MIS: 


"STe}04 Aep 6 - G wWoasT payertnoTeo suesay 
*pepunogz ere sazequmy 


eG, L°S 
ae 4°0 
ine rA Dai 
47° 0 S*O 
ak Chg 
cat z 
= care 
8°e L°4 
g°2 oe 
Or? T°S 
q M 


(°%,U09) OT°S aATEaVEL 


wis, 


It 


ueew 
YOOM OZ 


LY OT 
O€ 

Ez 

9T 
OT 


8 


_ ~ F¥* 


es 8 ae 


ae — 


203 


AIIM 

ONIGNS 
dOldsd JLVa 
- 5 
—— 4 
.g 
~----- M 
S 
= Sd 


02 
BOO 


él 


So AlVG NV3SW ‘81'S ‘Ol4 


o 


= - 


> 


er + Tr re ‘ 
esl? Set a ee, 


By July 22, the difference is reduced to a very small 
quantity. 

Although the above trend description may have to 
be qualified somewhat due to possible interpolation errors 
(see Section 5.3.5), there can be little doubt that F and 
G plots recorded higher Gs values than the remaining four. 
Figs. 4.6 and 4.7 reveal that from June to September 1975 
and 1976, F and G had the warmest 20 em and 100 cm soil 
temperatures, while the two wooded plots were the coldest. 
Preliminary indications point to Gs differences as an im- 
portant factor in the plots’ soil temperature differences. 
However, further analysis is required in order to quanti- 


tatively determine the significance of the Gs differences. 


DO. 3s H 

sensible heat flux data are listed in Table 517 
and illustrated in Fig. 5.19. H is the residual term in 
the energy balance calculations, since it was the only 
parameter that was not directly measured. Hence, its 
results are dependent upon measurements of Rn and LE, 
and to a lesser degree on Gs. 

Although there are large differences in mean values, 
with PS and S being the highest, and G the lowest, rankings 
for individual weeks varied. As with the LE results, trends 


in H are hard to see, but it appears that there is slightly 
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less chaos in the H curves so that certain patterns can 
be noted. In June, there is a general decline in values, 
followed by a sharp increase in July. A secondary minimum 
occurs on August 5, followed by another increase towards 
early September, and a decline of values at the end of 
the study period. However, certain individual plots do 
not appear to be following this trend at all times. G 
plot's values increase in mid-June, while the other plots 
showed a decline, and during late July and early August, 
B plot's curve slopes upward, contrary to the other five. 
The LE curve (Fig. 5.16) reveals that G plot's values for 
mid-June were fairly constant, so the increase in H was 
due to an increase in Rn alone. The situation in July 
represents a different story though, as LE values for B 
plot increased sharply, thereby accounting for the decline 
ary it. 

Other outstanding features include the large 
negative values for W at the end of June, and forG 
during the first week in August. Also, H values greatly 
fluctuate during mid-July and early August. The low values 
for W correspond to a period of very high LE values. The 
wide gap in H values on July 22 and August 12 were simi- 
larly caused. On the other hand, G plot registered only 


2 il 


a 30,000 Wm “day ~ increase in LE from July 29 - August 6, 


~2 -l 
while H values declined by more than 70,000 Wm “day ~. A 


decline in Rn accounted for the remaining 40,000 wn” “day~t, 
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In general, PS and S ranked high for most of the 
study period. W, F, and B exceeded the tree plots' values 
on several occasions, but were usually ranked between the 
tree plots and G, which had the lowest or near lowest 
values, except during late June and weekly periods ending 
August 12 and September 16. As with the other parameters, 
further analysis is required to determine the significance 


of these rankings before conclusions can be reached. 


5.3.4 Bowen Ratio 


Bowen Ratio calculations are listed in Table 5.12 
and illustrated in Fig. 5.20. The general trend during 
the first six weeks resembles the H curve, but instead 
of divergence, 8 values for all plots were quite similar 
at the end of June, despite large absolute differences in 
H. The tree plots separate themselves from the others 
during mid-July, while W does so on August 12 and F on 
August 19. There is incredible scatter of 8 values in 
September, but this is probably due to the decrease in 
magnitude of the energy balance components, which would 
facilitate enormous but unrealistic short-term changes 
in 8. Another possible explanation is the large tempera- 
ture difference between the warm surface and the fresh 
cold air mass that typically arrives in autumn. H (and 


consequently 8) would be greater than normal in this 
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Situation (Denmead, 1969). Subsequent calculations of 
mean 8 values should not be affected by this, in any case. 

In general, it can be noted that the early portion 
of the study period was dry, with more humid conditions pre- 
vailing at the end of June, and continuing towards early 
September. The exceptional values previously mentioned 
can be related to low LE values for these individual weeks. 
Low rainfall during much of July probably accounted for high 
8 values at PS and S, since light precipitation would have 
been almost completely intercepted by the tree canopies. 
More intense precipitation, such as that provided by storms 
experienced at the end of June and the beginning of August, 
are required before soil moisture levels benefit substantially. 
Since the first priority of energy use involves water use, 
more of the available energy could be utilized for trans- 
piration and evaporation rather than sensible heat. 

While the above may explain high 8 values for PS 
and S, the cause behind W plot's high 86 value in August 
remains a mystery. It is possible that there was an error 
in soil moisture calculations for this individual measure- 
ment. However, the detailed analysis that follows uses 
mean values compiled over the entire study (a period of 
20 weeks), which should effectively reduce the influence 
of this anomaly. Meanwhile, the high 8 value for F 


plot in August could be related to weeding and tillage, 
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which was done the previous week, and so temporarily 


reduced transpiration demands on Rn, thereby lowering 


LE, 


Ss eee, Errors and Preliminary Analysis 


Errors in energy balance calculations are basic- 
ally the same as those in radiation balance calculations, 
but with two major additions. The first involves errors 
in soil moisture calculations, and the second concerns 
interpolation of regression slopes for Rn and Gs. 

The neutron scattering equipment was calibrated 
during the study period, so large calibration errors are 
not anticipated. The most likely source for measurement 
error was the placement of the surface probe. This created 
some difficulty, especially at W plot, because the shrubs 
themselves were fairly close together. Since the August 
12 measurement at W plot cannot be easily explained by 
physical factors, it may be assumed that the surface probe 
was not positioned properly on the plot surface. APW 
calculations were based on soil moisture totals which 
were weighted according to the measurement depth. For 
example, if data were available for a 1 metre depth, the 
calculation procedure would be as outlined in Table 5.13. 
The surface calculation thus has a 10% influence on the 


final sum if data are available for 1 metre, and a some- 
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TABLE 5.13 
PROCEDURE FOR SOIL MOISTURE CALCULATIONS 


Listeall datasby depths» e.g00e0s), mee-exl 


0.2 m -- x2 
0.5m -- x3 
0.8 m -- x4 
1.0 m -- x5 


Calculate mean values for each layer; e.g. xl + x2 


Weight each value according to width of soil layer; 


See Re oe kee 53 x 0.3 


Sum all values to produce weighted mean; units are 


fm soil 
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-l -1 s nt 
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what greater influence if only shallower depths can be 
monitored (as was the case with PS and S plots). This 
appears to be an isolated case, however, as the various 
soil moisture curves (including Fig. 5.15) show W following 
trends similar to several other plots both before and after 
this*date. Final analysis is based on 20 week mean values, 
and it is anticipated that an individual error such as this 
one will not significantly affect the final results. 
Regression equations have already been discussed 
in this chapter. Several large time gaps exist between 
periods of direct measurement, so there may be errors in 
interpolation of regression slopes during these periods. 
Due to the high correlations and significance levels re- 
corded for all regression equations used, the magnitude 
of this error should be small, although this would be dif- 
ficult to verify. Similarly, the usage of eq. (2.6) and 
assumed values for Gs/Rn poses additional problems on 
late August and early September results. However, as 
with the aforementioned situation, the regression equations 
utilized show high correlations and significance levels, 
and since Gs is a very small quantity, absolute errors 
should be small for both estimates. As previously men- 
tioned, the final analysis is based on 20 week mean values. 
Also, shorter period mean values are calculated for 2 week 
periods to facilitate the equalization of time intervals 


for all measurements at all plots during the entire study 
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period. In other words, where there had previously been 
periods of 5 - 9 days, the final analysis will use only 

14% day intervals. This should also serve to reduce the 

effect of random fluctuation errors, 

Energy balance percentages for each period are 
‘Shown in Table 5.14, Dates where data were estimated 
are included. Besides those previously mentioned, there 
were days when individual hours were missing from the base 
station data, but these were few in number, and should not 
affect mean values. The assumed Gs/Rn values are very 
small in magnitude, except for F plot, but large errors 
are not expected for this case either, since 15 weeks of 
data are available upon which Gs/Rn estimates can be made 
with some confidence. 

[he data presented in this study have been shown 
to be reasonably accurate, and representative of actual 
field conditions. Having considered the various error 
possibilities, it can be concluded that the data can be 
analyzed at least on a comparative basis. The purpose 
of the statistical analysis to follow is to determine 
the significance of the differences in energy balance 
between the plots. Upon examining mean energy balance 
ratios (Table 5.15), differences as large as 0.2 are 
noted for H/Rn, and 0.5 for 8. It is hoped that when 
these differences are examined, the picture given by 


this table will be made clearer. 
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5.4 Statistical Analysis of Energy Balance 
5.4.1 Analysis Procedure 


The purpose of this analysis is to determine trends 
and evaluate the statistical significance of the differences 
between mean values of the energy balance components for 
each plot. Having established that the data reasonably 
represent field}conditions, the Paired * test ican! be 
used to provide an index, t, which is the ratio of the 
difference between the means to the standard error of the 
difference. This procedure examines sample means to see 
if the differences between them were a result of random 
variation, or significant differences between the individ- 
ual cases that constitute the mean. This assumes that the 
individual differences are normally distributed. Several 
histograms illustrating the distribution of these differ- 
ences are illustrated in Fig. 5.22%. In a normal distri- 
bution, 95% of the individual cases are found within 2 
standard deviations (SD) of the mean. This condition 
is satisfied in 7 of the 8 distributions shown. When 
dealing with small samples of data, certain inaccuracies 
can result from skewed distributions, and these must be 
taken into account. 

The individual cases are 2 week totals of Rn, H, 
LE, and Gs, for each plot. This was chosen to equalize 


the cases' time intervals. The t test is a rigorous 
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FIG. 5.21: HISTOGRAMS 
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technique, and can be used for sample sizes as small as 

3 cases (N = 3), if the concept of "degrees of freedom" 
is employed, rather than the actual number of cases. 
Degrees of freedom are the number of values that can be 
assigned arbitrarily, assuming that the sample total and 
mean remain unchanged. Although N-1 is often used for 
determining degrees of freedom, N-2 is used here since 
the standard error is also assumed to be unchanged. This 
increases the difficulty in proving significance, since a 
higher t value must be attained before the critical value 
is reached. For the purpose of this analysis, the 95% 
confidence limit is chosen as the minimum for rejection of 
the null hypothesis (there is no significant difference 
between the samples) and acceptance of the alternative 
hypothesis (there is a difference). 

Simultaneously, the data are analyzed for corre- 
lation (r). This test determines the degree of relation- 
ship between the trends of a pair of samples. High cor- 
relations imply a strong statistical relationship, with 
a perfect correlation being r= 1.0. Note that it is 
possible to have high correlations and insignificant t 
values, and vice versa, for the same pair of samples, 
since the tests are independent of each other. 

Before examining the test results, it is important 
to point out the limits of what a statistical test can say. 


Although it is possible to prove that the means are signifi- 
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cantly different, or that a strong correlation exists 
between two samples, the interpretation of these results 
is essentially based on human judgement and experience. 
These tests provide a mathematical basis for decision 
making, but they are not an end in themselves. Human 
judgement requires examination of many factors besides 
the test results, including physical, climatological, and 


biological factors, before conclusions can be reached. 


5.4.2 Results 


Paired t test results for Rn, H, IE, and Gs, 
are listed in Table 5.16. Analysis was done for the 
entire 20 week period, and for the first and second 10 
week periods separately. Those results significant at 
the 99% confidence limit, as well as those at the 95% 
and below 95% (not significant) levels are included. 

The numbers in the matrices are differences between 
the means of the two samples indicated. The difference is 
calculated by subtracting the row sample from the column 
sample. Negative numbers indicate that the row sample 
had the higher mean value of the two. 

The rows and columns of each matrix are ranked 
by the magnitude of their sample mean daily Rn values. 
This parameter was selected because it is dependent on 


all other energy and radiation balance factors. It was 
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found that PS had the highest value, followed by W, F, S, 
B, and G, in that order. Since this ranking is used in 
every matrix, the negative values indicate that the 
ranking of that parameter for that time period does not 
agree with the 20 week (or total period) mean daily Rn 


rankings. 


Correlation test results are listed in Table 5.17. 


The arrangement of these matrices is identical to Table 
5.16, as are the indicated confidence limits. As pre- 
viously mentioned, the distribution of significant t and 
r relationships are not necessarily the same. 

At first glance, it can be noted that for the 
total period, significant differences exist between the 
plots' Rn and Gs values. In addition, strong degrees of 
relationship exist between the plots. However, if the 
period is split into halves, the first half shows poor 
correlations for Gs, while most of the second half dif- 
ferences are found to be not significant. Upon examining 
H, it appears that G plot is significantly different from 
the rest during the total and second half periods, while 
certain other plots, notably W and F, show significant 
differences from the trees and the remaining open plots 
respectively. Differences in the plots’ LE values are 
found to be not significant over the total period, while 


some differences are computed between W and the trees 


during the first period, and B and several plots during 
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TABLE 5.16 
T TEST RESULTS (1000 Wm7day~+) 
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TABLE 5.16 (con't) 
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TABLE 5.16 (con*t. > 
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TABLE 5.16 (con't. ) 
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the second. Both H and LE values show large numbers of 
significant correlations over the total period, but there 
are almost none for H during the second half. 

Another general observation that can be made 
concerns the distribution of negative differences. Rn 
and H have few, but LE and Gs have many, indicating that 
their rankings are considerably different from the total 
period mean daily Rn rankings. Having based the plot 
rankings on total period mean daily Rn values, it is ap- 
propriate to concentrate the analysis on this particular 
matrix. The other matrices can provide additional input 
as to the factors behind the establishment of this particu- 
lar ranking, and the implications of these results on micro- 
climate. 

According to the test results, PS has significantly 
higher Rn values than any other plot. There appears to be 
little difference between W, F, and 5S, so they may be 
grouped between PS and B. G plot is significantly lower 
Thanes in the Kn ranks: 

Why is the PS plot ranked first? It was found to 
have relatively low albedo values, but S and F plots' al- 
bedos were comparable. However, PS had lower longwave 
losses than S and F, higher first period H values than 
F, and higher total period H and Gs values than sS. Its 
second period Rn values were considerably higher than F 


and S plots’ values. F was not significantly lower during 
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the first period since its low weed cover allowed greater 
soil exposure, thereby reducing albedo. Meanwhile, what 

of W? This plot's higher albedos were compensated by lower 
longwave emissions, while its low September albedo allowed 
greater radiation absorption during the second period. Non- 
significant differences were found for the two half periods, 
yet the situation is different for the entire period as a 
whole. The only source for this difference is first period 


H, where PS is higher by 54,000 wim“ - 


day This is partially 
offset by significantly higher values for first period LE at 
W plot. Both relationships have strong correlations, as do 
the Rn differences. A re-examination of first half field 
data reveals that W had much greater soil moisture content, 
and LE demands were higher than the Rn available, thereby 
necessitating additional input of advected energy from 
elsewhere. Consequently, H became negative. Had W plot's 
albedo been lower, Rn might have been raised sufficiently 
to meet LE demands and so reduce the Rn difference. Hence, 
the key factors behind PS outranking W are probably lower 
soil moisture content and lower albedo during periods of 
high LE demand. 

It is important to note that PS plot's LE values 
@id not differ significantly from FP, and during the total 
period, a fairly strong correlation was computed. This 
indicates that despite differences in age, size, and soil 


moisture availability, PS and P plots' seasonal evapotrans- 


piration rates are similar. However, when P was compared 
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to W, the other deciduous plot, it was found that P plot's 
IK values. were significantly higher .than W plot's during the 
second period. .P is older and larver than the other two, 
and this should reduce its ability to draw in heat 

and energy from the outside, in which case its LE values 
would be reduced somewhat. Differences in species, age, 
plot size,, and soil moisture availability, are all factors 
in LE differences. However, in this case, these factors 
have combined to reduce the LE differences between PS, P, 
and W over the total period. Otherwise, it is likely 

that PS, with adequate soil moisture, would have far ex- 
ceeded the other plots' LE values, and its H values would 
have been reduced. It is doubtful that the Rn rankings 
would be significantly changed, however, because of the 
radiation balance differences, especially the albedo 
component. 

The next question concerns W, F, and S plots’ 
similar Rn rankings. A significant difference was found 
between W and S during the total period, and W and F 
during the second half. No significant differences were 
computed for the first half. Correlations were highly 


Sicnificant sin all wases. 


lonis is the "oasis" effect (see Rauner, 1976), which 
was previously discussed in Section 2.4.1. The effect 
becomes weaker as woodlot sizes increase. 
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Why should these plots have fairly similar 
rankings? W plot's albedo is much higher than the other 
two while it has a leaf canopy. However, longwave losses 
are lower, so W plot's Rn values are raised. The source 
of the difference between W and S appears to be total period 
Gs, but nothing seems to explain the 9,000 Wm7 “day advant- 
age calculated for W over F during the second period, since 
all other components show non-significant differences and 
poorly correlated relationships. The total period result 
shows no significant difference between W and F, so the 
second period's difference was not large enough to affect 
it. There is probably a trade-off between F plot's advant- 
age in Gs, and smaller less significant advantages in H and 
LE for W and S. It would appear that W plot's higher albedo 
is compensated by greater evapotranspiration rates and lower 
longwave emissions. The latter is probably a result of the 
former, since W plot could have had a higher relative humid- 
ity on a micro-scale, which would tend to reduce longwave 
losses. On the other hand, trees are known to reflect long- 
wave radiation more strongly than solid ground (Geiger, 1965; 
Jarvis et al., 1976). Therefore, differences in canopy 
structure, bare soil exposure, and foliage characteristics 
(e.g. leaf orientation) must be taken into consideration 
when examining Ln losses, as well as for the Gs advantage 
enjoyed by W because of its stand structure, which is lower 


and less dense than S. 
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The two lowest ranks are occupied by B and G plots. 
The magnitudes of the Rn differences between B vs G and B 


2 aay t respectively. 


vs S are 20,000 Wm” “day7t and 9,000 Wm 
Both values were found to be significant, and the relation- 
ships well correlated. 

B plot undergoes rapid structural changes during 
the growing season, as it advances from seedling to mature 
plant in four months or less. Consequently, two important 
parameters are affected; albedo, and soil water use. Al- 
bedo varies with plant height and bare ground exposure 
(Fritschen, 1966), so reflectivity was higher for July's 
Tr rleceplany COVer, tian LOn uune Ss GhOrvecrmOne. ns wie 
plants matured and changed colour, reflectivity was reduced 
substantially. LE values were high from June to August, 
but declined markedly in early September, due to reduced 
plant activity. In fact, second period LE was shown to 
be significantly lower than G plot's, while its H values 
were higher. Second period H appears to be one main factor 
behind B plot's higher rank, since G plot had higher total 
period Gs, as well as second period LE. Another factor 
that must be considered is longwave loss, which is greater 
for G plot fram June to August. 

Meanwhile, what about the plots ranked above B? 

It appears that albedo is the main factor separating B 
from W, F, and S. The difference between B and F is aug- 


mented by F plot's greater Gs values. W had lower longwave 
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losses than B, especially in June, perhaps indicating 
lower micro-scale relative humidity and higher surface 
temperature at B. The smaller Rn difference between B 
and S is due to albedo differences alone, and this is 
slightly offset by a.2,000 wm72aay~+ advantaze for.B plots 
Gs values. Previous studies in Canberra, Australia (Denmead, 
1969) indicated that Rn and LE should be higher for conifer- 
our species, while H and Gs are lower. The discrepency in 
the Ellerslie study's H and LE results must therefore be 
caused by 58 plot's limiting soil moisture supply. 

G plot was ranked sixth, due to its high albedo 
and longwave losses, and significantly lower H values. 
These were slightly offset by high Gs values. Lé& values 


did not differ significantly from other plots. 


SF Implications of Results on Central Alberta Climate 


eer en i eS ca en 


ae Basis for Conclusions 


Previous sections have included discussions on raw 
field data, calculated energy balances, computed statisti- 
cal results and the reasons behind them. Various errors 
have been considered and interpretations presented within 
the limits set by the quantity and accuracy of the data 
and the power of regressions, t tests, and correlations. 
As a result of these factors, there are no attempts at 


utilizing absolute values of energy balance parameters 
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for the purpose of mathematical or statistical modelling 
Similar to Table 2.3 or otherwise. However, there exists 

a statistical basis for comparing the gross energy balances 
of the various plots at Ellerslie. That was the purpose of 
the field program. It is now possible to discuss the impli- 
cations of these results on Central Alberta climate, keeping 


in mind the potential influence of horizontal advection. 


5 one Potential Micro-Scale Effects 


Significant differences in Rn and Gs have been 
calculated for all six experimental plots. Certain plots 
also exhibited differences in H. No significant differences 
were calculated for total period Li although some were found 
for shorter periods. Since all field investigations were 
done on relatively flat terrain, any statements made here- 
in will apply strictly to flat areas. Large topographic 
variations have a profound influence on energy and water 
balances but exposure and slope aspect were constants in 
this experiment, and are not considered in the following 
discussion. 

PS has the highest Rn values, so one would expect 
that its surface temperature would be relatively high as 
well, yet its longwave emissions were the lowest of all, 
despite what should have been a high longwave albedo. 


Deciduous trees are efficient water users, so even though 
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soil moisture levels were low, the trees survived the summer 
with no apparent damage. However, its high H and g values 
are interesting, especially since LE siteaa normally be 
dominating its energy balance (Rauner, 1976). Since H 
increases if surface temperatures are warmer than the air, 
it@is’ possible’ that°H-is most =significantvat night, wien 
air temperature and evaporation rates are lower. Conse- 
quently, a land use change from farmland to deciduous 
forest would probably cause an increase in nighttime mini- 
mum temperatures at canopy level. If soil moisture is not 
limiting, daytime LE should increase to the point where it 
approaches, or maybe exceeds Rn. In that case, daytime 
maximum temperatures would be reduced while relative humid- 
ity would rise. This effect would be increased if there 
were numerous woodlots of small size (less than 1 xm“ ) in 
fairly close proximity to each other. Unfortunately, soil 
moisture is frequently limiting, so chances of this occurring 
in Central Alberta are low, unless several successive wet 
years occur. In either case, soil temperatures will remain 
lower than the more open surfaces due to significantly lower 
Gs values. 

Coniferous woodlots should have slightly cooler 
soil temperatures than deciduous areas, since Gs is lower. 


However, other conclusions involving coniferous areas are 


not as easy to make. 
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Longwave losses are higher, implying either a higher 
longwave albedo or higher surface temperatures. H and Rn 
are significantly lower for the total period and these re- 
lationships are highly correlated so the trend is consist- 
ent. Since lower H and Rn values imply cooler surface temp- 
eratures relative to the overlying air, longwave albedo 
would ordinarily be the most likely cause of high In losses. 
Although energy balance results would seem to indicate that 
coniferous areas should be cooler than deciduous areas, 
screen level temperatures were warmer at S plot due to 
canopy structural differences. On the other hand, lower 
H values may be due to an increase in overlying air tempera- 
ture, rather than a decrease in surface temperature. The 
situation is further complicated by its extremely low soil 
moisture levels, which affected LS. 

There is one factor that has not yet been discussed; 
Rn variations with surface wetness. Wilson's (1976) study 
of irrigation effects on micro-scale climatic parameters 
at Ellerslie showed that Rn was lower over dry grass than 
wet grass. He attributed this to an increase in surface 
temperature, which would increase longwave emissions there- 
by reducing Rn. Therefore, it seems likely that 5S plot has 
higher surface temperatures than PS, so Vt sLtSesoiverconal— 
tions continue to be drier, coniferous area air temperatures 
will be warmer than deciduous. Since H values are lower, 


these must be attributed to lower Rn, and in an indirect 
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way, to low soil moisture levels. g does not differ sig- 
nificantly from PS plot's results (Table 5.14) so small 
coniferous woodlots will have an effect similar to small 
deciduous woodlots if soil moisture is enhanced. As long 
as the prevailing weather remains dry, its “oasis” influence 
will be reduced. Unlike the deciduous areas, coniferous 
areas would probably exhibit higher daytime and lower night- 
time air temperatures, anda lower relative humidity than 
both deciduous and open areas. It would be interesting to 
compare S and PS under identical soil moisture conditions. 
The shrub plot ranked second in mean daily Rn, 
primarily due to low longwave losses. Although it is a 
deciduous cover, its unique canopy structure is responsible 
for creating its own micro-climate, which is different from 
a deciduous forest. Soil moisture levels were much higher 
at W than at PS, because of its superior ability to retain 
winter snows. This additional snow cover was instrumental 
in maintaining higher winter soil temperatures. However, 
the high moisture content probably raised the soil's thermal 
conductivity, since water replaced air in the soil pore 
spaces. In addition, there was less organic matter than 
at the other open ete Consequently, W had lower summer 
soil temperatures than B and G plots, despite a lower albedo 


and comparable Gs values. Therefore, shrubs should maintain 
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warmer winter and cooler summer soil temperatures than 
fallow, grass, and barley covers. Its soil moisture levels 
should continue to be higher as well, although its advantage 
could be reduced after several dry winters and wet summers. 

It is somewhat difficult to state any firm opinions 
about air temperature over shrub covers. No direct air 
temperature measurements were made, as these were assumed 
to be similar to G plot's. Due to the small size of the 
plots, this assumption was valid for the experimental con- 
ditions, but what about for larger surface areas? Results 
indicate that Ln is lowest and Rn is highly ranked, but 
total period H does not significantly differ from any plot. 
Since its soil moisture levels were consistently higher, 
its surface temperature was probably cooler than Ps, 8, 
and F plots, while a full leaf canopy existed. Lower al- 
bedo, similar or higher In losses, and higher H values, 
indicate that autumn surface temperatures may be as high as 
or higher than these plots. 

g was higher at W plot than at B and G plots, but 
statistical analysis reveals that total period H and L& 
differences are not significant. Therefore, W plot's sur- 
face temperatures are probably cooler than these plots as 
well, when it has a full leaf canopy. During autumn, the 
situation is similar to that previously described. These 
statements raise questions about differences in time of 


leaf fall between PS and W, and the changes in late summer 
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and early autumn energy balances that this causes. 

The present soil moisture regime is the main factor 
behind these differences. Higher albedo contributes to the 
difference between W and PS, S, and F plots during fuil 
canopy conditions. Since shrub's albedo is lower than B 
and G plots' throughout the growing season, soil moisture's 
influence increases in this situation. Another factor that 
may be considered in the latter cases is lower wind speeds 
over shrub than over barley or grass. This’is’a likely 
event, considering shrub's taller woody canopy structure, 
which would offer stiffer resistance to surface winds. The 
jdea is speculative, and a comparison of screen level wind 
speeds and roughness length calculations would be useful in 
determining its accuracy. 

F plot ranked very close to W plot in Rn, but its 
radiation and energy balances were quite different. Higher 
longwave losses were evident despite lower albedo values 
during the first period. Gs was much higher than at all 
other plots due to greater exposure of the dark soil sur- 
face to incoming radiant energy, so higher soil temperatures 
should be expected under fallow fields during summer. 
Winter temperatures will vary with depth of snow cover, 
and under existing conditions, lower temperatures are ex- 
pected unless there is a heavy weed cover or inducement of 


snow drifting by human intervention (e.g. snow fences). 
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Despite its advantages in Rn, F plot's H values 
were not significantly higher than barley's, except during 
the first period when barley was actively transpiring. G 
plot's H values were lower for the total period, while no 
significant difference was calculated between F and PS, S$ 
or W. 

Having concluded that F plot's summer soil tempera- 
tures are highest, it must now be determined whether this, 
as a surface temperature, is higher than the other plots' 
surface temperatures. It appears that the answer is yes, 
since PS, W, and B plots have lower Ln emissions. However, 
this is not obvious when comparing F to S andG plots. 

S plot recorded comparable longwave losses because 
Rn and albedo were not very different, and on occasion, its 
losses were actually greater. Overall, their surface temp- 
eratures are probably similar, although this relationship 
is certain to be affected by weed cover density and surface 
wetness at fallow fields. 

Similar longwave losses were also calculated for G 
plot, but Rn is much lower, and albedo much higher. Mean- 
while, soil temperatures under grass were higher from July 
to September despite lower Gs values (see Fig. 4.6). G 
plot had greater water use, but F plot had lower soil mois- 
ture levels at the start of the study period, which may 
have reduced LE differences. It is suggested that F has a 


higher daytime surface temperature than G because of lower 
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albedo and soil moisture levels (which ultimately lead to 
lower evaporation rates). At night, these higher moisture 
levels, combined with the overlying plant cover, should 
allow G to maintain warmer soil temperatures, but the air 
layer in contact with the surface would still be cooler 
than at F plot since release of heat is more gradual. As 
a result of higher daytime evaporation and transpiration, 
relative humidity of air over grass is greater, and when 
combined with cooler nighttime surface temperatures in 
September, dew occurred more frequently over grass. Since 
maximum dew usually occurs 1 - 2 metres above ground in 
non-forest areas (Geiger, 1965), the taller grass plants 
would have an added advantage in attracting dew. 

Thus, the overlying air is cooler over grass than 
over fallow fields, and this temperature difference should 
increase if the fields are larger in size (as opposed to 
the reduced “oasis" effect should woodlot sizes increase). 
Relative humidity will be lower for fallow field areas, 
except during immediate post-rainfall periods, in which 
case surface evaporation would be higher and 8 would de- 
crease. Subsequently, air temperatures would be lower, 
and relative humidity would be higher. 

One question that remains is whether G plot's 
surface temperature is cooler than all other plots besides 


F and S. PS and W have previously been examined. It was 


concluded that PS would have higher nighttime minimum 
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temperatures, and possibly lower daytime temperatures, 
depending on LE. W would be cooler when its canopy is 
full, but might be warmer after leaf loss occurs. 

B plot had higher Rn, lower Ln, and lower Gs values 
than grass. Second period H was much higher while second 
period LE was lower. Overall g was higher, but this result 
was still below Ps, S, W, and F plots. Since LE appeared 
to dominate B plot's energy balance from June to early 
August, its surface temperatures were probably cooler. 
second period results indicate reduced water use, so more 
energy was utilized as sensible heat. Therefore, barley 
surface temperatures during second period were probably 
higher than grass temperatures, and the overlying air's 
relative humidity was probably lower, which reduced the 
frequency of dew formation. However, it would not exceed 
temperatures at the fallow and tree plots. This does not 
agree with Denmead's (1969) findings, which showed (Pinus 
radiata) tohave lower surface temperatures (mean for the 
canopy) than wheat, since Li was higher. At #llerslie, 
however, soil moisture is extremely limiting at 5S plot, 


so LE is reduced and H increases. 


Sey se SI Potential Macro-Scale #ffects 


Rural land use in Central Alberta is dominated by 


cropland. Census figures indicate that this dominance is 


anuiev gp ‘awed - Sia en 
broveq ef dstw, cadlginl doue sew W bolasd baade®’/eeeny hat 
siueer abte ied ytedgla’ caw goiferev® sctawol'ssw BD boiteg 
basgoqqea Sy eonle: petolg F bre yw S ,e&:wolod [Lites esw 
yliee of shphi mort: wotealsd yerete: wi voly & otaminod oF 
aslooo. gidsdez¢q oxow; gevutbiagmet eontse: etl <\teumA 

stop oe ,@ap tevaw heovbex atesibal-atlses: bodreq baoose 
vetyed ,exobevad? «,taed sidieres. 26 bosility tew vetete 
cuiteden quoncbelueq! basen! soisuhi aetussonqmat poetIve 
e'sis poiyiteve edf- bas ,eoratarsqmsy Sesta- -padid -redgit 
ont hoguber toitw.4 tewol Yi cedate Baw: ed lbamur evitales 
beepka) tom biuoy $2 vtevewol =, ft £SnarsoT ‘woh ‘re° yonenpert 
ton noob aint. .cetedg sow bre. wodlat: cht ta) eatute croquet 
evel?) pewode doddw \emrtbatthe(eoenye etheoandd dtiw-sotge 
ett 20% ser) 2eqteraqmet enetiie towel etad oF (atsibs 

.  otierella sa. cretgid:eaw Gt eorie: ;teedwinedd (yqonse 
aolq Bote gaetinis Ylomertxeval \ewtetemtlos sevewor 
peseeetosi KH bia beoubet oi) alos 


a 2 
ie 


os havantoon wh nénesia tasted al ous bank tari pare 
th sonatinyd abit dad ofsoiinl somugit eveded, -basiqexe 


seo 


; 
; 
: 


increasing, as is the use of land for fallow rotations. 
Natural pasture and woodlots are declining. It is beyond 
the scope of this discussion to speculate on future rural 
land use trends. Some hypothetical possibilities can be 
mentioned, however, and their associated effects on climate 
described. 

If woodlot acreages continue to decline in favour 
of cropland and fallow fields, albedo levels in former 
wooded areas will show an increase in general. Surface 
level temperatures should increase, while at former canopy 
(tree-top) levels, temperatures may decline slightly. Soil 
moisture shall be enhanced since interception of winter 
snows will be reduced, but surface wind will increase, 
and unless shelterbelts are present, much of this snow will 
be blown away. There will be other horizontal effects as 
well, including an increase in open area horizontal heat 
transport. Surface temperatures will depend on the extent 
of fallow acreage, but winter and summer soil temperatures 
will certainly increase. If soil moisture conditions are 
improved, LE will Fe EPS erie a continuation of present 
conditions will lead to increased H values, higher surface 
and air temperatures, and wind erosion will become more 
prominent. Vertical effects include increased temperature 


gradients, which could lead to greater local convective 
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On the other hand, should woodlot acreages increase, 
surface winds will be reduced, albedo will decline, and soil 
moisture and soil temperature will decline. 8 will increase in 
woodlot areas, but their horizontal "oasis" effects could 
be profound, if, as previously described, numerous small 
woodlots are established in the region. However, this 
depends on availability of adequate soil moisture and it 
is anticipated that under present conditions, soil moisture 
will not be adequate under flatland woodlots.~ Shelterbelts, 
however, are a different matter. Rows of trees oriented 
perpendicular to the wind can act as snow fences so that 
Snow can accumulate. Open shrub woodlots would also be 


useful in this regard. 


Ipesides interception, woodlots may prevent the accumu- 
lation of snow at the soil surface by enhancing snowmelt 
and evaporation. This would be a consequence of high 
radiation absorption at the woodlot's edge (Hage, per- 
sonal communication). 
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CHAPTER VI 


SUMMARY AND CONCLUSIONS 


6.1 Summary of Findings 


The gross energy balances of six different land- 
scape surfaces typical of Central Alberta have been ex- 
amined. Analysis included measurements and derivations 
of albedo, net radiation, sensible heat flux, latent heat 
flux, soil heat flux, and other parameters such as soil 
temperature, soil moisture, and air temperature. Incoming 
global solar radiation, precipitation, and topographic 
exposure were constant for all plots. 

Results have been given for each of these parameters, 
and comparisons made between the different plot surfaces. 
The plots have been ranked according to total period mean 
daily net radiation. The small poplar plot ranked highest, 
followed by shrub, fallow, spruce, barley and grass DLOtS 
in that order. The shrub, fallow, and spruce plots were 
found to have fairly similar net radiation values, while 
significant differences were calculated between the re- 


maining surfaces. Therefore, a more appropriate ranking 


would be: 
highestes..) = DOpDLar 
shrub, fallow, spruce 
barley 
lowest .. . §rass 
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The major factors behind this ranking appear to be dif- 
ferences in albedo, canopy structure, relative humidity, 
soil moisture, and longwave emissions, but these had 
varying influences on each of the net radiation results. 
No significant differences have been calculated for total 
period latent heat flux, indicating similarities in total 
period water use brought about by limitations in soil 
moisture. Grass had the lowest sensible heat flux values 
over the total period. Fallow was found to have the highest 
soil heat flux values, with grass ranked second. Upon 
partitioning the study period into halves, significant 
sensible heat flux and latent heat flux differences were 
calculated between several plot pairs. Therefore, it is 
not surprising to find that hypothetical rankings of sur- 


face temperatures would probably vary with time: 


midsummer late summer, early autumn 
highest... fallow, highest... fallow, 
spruce spruce, 
shrub 
poplar poplar 
barley barley 
grass lowest... grass 
lowest «6 eee enrud 


Experimental results are not in complete agreement 
with published findings from other locations. Areas of 


discrepancies include energy balance results for the two 
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wooded plots. Specifically, latent heat flux values are 
below expected levels, while sensible heat flux values 
are above. It is concluded that limiting soil moisture 


is the main factor behind these changes in energy balance. 


(ie Conclusions and Recommendations 


The purpose of this study was to examine the gross 
energy balances of six different vegetation surfaces, and 
the manner in which each surface responded to equal amounts 
of incoming solar energy and precipitation. Despite the 
small size of each plot, significant differences were iden- 
tified within their respective energy balances, and these 
were linked to several unique factors found in the various 
vegetation covers. 


Net radiation is not a purely atmospheric factor. 


Its magnitude may be strongly affected by surface properties 


such as albedo, canopy structure, and soil moisture avail- 
ability. Evapotranspiration is similarly dependent on 
canopy characteristics, as well as macro-scale atmospheric 
events. Consequently, both factors have to be considered 
before applying published energy balance data from other 
locations to local conditions. 

Differences between external published energy balance 
data and local data are also indicative of the adaptability 


of certain vegetation covers to environments that may be 
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outside their normal ranges (or habitats). Ellerslie is 
located in an aspen transition zone bordering natural 

prairie grasslands. The two tree plots, especially the 
coniferous, have unusually low soil mositure contents and 
high Bowen Ratios. They are absorbing too much energy, and 
this is probably affecting their surface temperatures. Nor- 
mally, local deciduous and coniferous woodlands are restricted 
to stream valleys and depressions where there is some shade, 
and moisture can accumulate. Large scale tree planting on 
open flat lands will not be very successful unless additional 
moisture is provided, and tree densities are restricted to 
fairly low levels. The coniferous plot at Ellerslie is too 
dense for the local environment to adequately support it. 

It is likely that the spruce trees here will not survive 
much longer, unless several wet summers occur. Individual 
trees and shelterbelts are more suitable than square wood- 
lots. 

Many questions arise out of this study. If net 
radiation differences can be traced to differences in sur- 
face temperature between wet and dry surfaces, does dur- 
ation of surface wetness vary between the plots: is the 
deciduous canopy warmer or colder than the coniferous can- 
opy? Would their albedos and air temperature profiles be 
similar if soil moisture is unlimiting? Would greater soil 
moisture levels affect their soil temperatures? What are 


the radiation and wind profiles of the tree and shrub plots? 
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Do these plots affect micro-scale wind patterns? Would 
changes in wind patterns and the differences in canopy 
structures influence snow accumulation? Are there sig- 
nificant “edge effects" at the tree plots which could be 
partly responsible for reducing their soil moisture levels? 
Would their effects be different if the plot sizes were 
larger? Does tilled fallow use more water than weed in- 
fested fallow? Are the screen level air temperatures at 
the shrub, barley, and fallow plots different from the 
grass plot? Would they be different if the plots were 
larger, and if so, would this affect evaporation, and 
consequently, plant growth and soil moisture levels? 

Would energy balance results be significantly different 

at other locations within the area, such as Stony Plain, 
St, Albert, or Wetaskiwin? Research on these and other 
topics are required before land use influences can be more 
adequately predicted. 

This study, despite its size, has merely scratched 
the surface. The Ellerslie facility is unique, and should 
be utilized to its fullest capacity. In fact, it should 
be expanded if possible, to include plots of different 
crop types, and additional tree plots, since mean values 
could then be calculated from several plots rather than 
just one. This would improve data. accuracy immensely so 
that absolute values of net radiation, latent heat flux, 


soil heat flux, and other parameters could be utilized 
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directly for the purpose of developing a statistical 
model that could predict changes in energy balance from 
changes in land use. Facilities of this kind should also 
be established at other locations if possible. 

Geiger (1965) quotes the late German researcher 
A. Schmauss, who said that “everything fits in with 
everything else." This is particularly true when the 
living processes of plants are involved, since there is 
always feedback between plants and environment, and this 
feedback may affect the physical relationship that is 
being examined. The adaptive ability of plants must be 
considered when attempting to assess land use effects on 
climate. But plants have limits to their ranges, and 
human intervention is often necessary to maintain these 
plants outside their natural habitats. The present dilemma 
that many North Americans are facing this year is partially 
a result of inadequate consideration of local climatological 
and hydrological resources, and the limits they impose on 
local agriculture. No extensive alteration of the land- 
scape should ever be undertaken without first accepting 
the responsibility for any side effects this alteration 
may cause, whether or not a problem might arise "only" 
once in twenty years. Since these side effects are not 
yet adequately documented, for Central Alberta at least, 
such action should not be blindly undertaken or should be 


adequately documented when action has already been undertaken. 
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